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METHOD AND APPARATUS FOR DETERMINING THE HEALTH OF A 
COMPONENT USING CONDITION INDICATORS 



Background of The fnvenfinn 

1- Field of the Invention 

This application relates to the field of vibration analysis and more particularly 4 
pafoiming vibration analysis for the puipose of device monitoring. 

2. DescriPtioii of Relatftd Art 

The transmission of power to rotors which prx)pel helicoptere and other shafts that 
propel devices within the aircraft induce vibrations in the supporting structure. The 
vibrations occur at frequencies that correspond to the shaft rotation rate, mesh rate, bearing 
passing firequency, and harmonics thereof. The vibration is associated with transmission 
error (TE). Increased levels ofTE are associated with transmission faihire. Similar types of 
vibrations are produced by transmissions in fixed instaUations as well. 

Parts, such as those that may be included in a helicopter transmission, may be 
replaced in accordance with a pr^ietennined maintenance and parts replacement schedule. 
These schedules provide for replacement of parts prior to failure. Tht replacement schedules 
may indicate replacement time intervals that are too aggressive resulting in needless 
replacement of working parts. This may result in incurring mmecessary costs as airplane 
parts are expensive. Additionally, new equipment may have instaUed faulty or defective 
parts that may fail prematurely. 
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Thus it may be desirable to provide for an efBcient technique for detectingpart and 
devicedegradationwithoutunnecessarilyreplacingparts. It may be desirable that this 
technique also provide for problem detennination and detection prior to failure. 
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In accordance with one aspect of the invention are a method executed in a computer 
system and a computer program product for determining a health indicator associated with a 
5 component. A phffaUty of health classifications are detennined. At least one condition 
indicator is detennined quanting a characteristic of the component. A probabiUty 
associated with each of the health classifications is determined. The probability is an 

estimation that the component is of a particular health classification given the at least one 
mdicator. A detennination is ftiade as to which of said health classifications is associated 
10 with said component using said probabilities associated with said health classifications for a 
given set of observed values. 

In accordance with another aspect of the invention are a mefliod executed in a 
computer system and a computer program product for deiermining a health status of a 

15 component A plurality of condition indicators are selected having a vahie and each having a 
corresponding weighting factor, and at least one threshold value defining at least two 
classifications. A contribution to a health indicator is determined for each of said condition 
indicators, wherein said determining fiirther comprises, for each of said plurality of 
indicators: determining which of said at least two classifications said value of said each 

20 indicator belongs; and determining said contribution to said health indicator by said each 
condition indicator in accordance with a selected one of said at least two classifications and 
said weighted value. The health indicator is determined in accordance with all contributions 
by each of said condition indicator values. 



3 
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In accordance with one aspect of the invention are a method executed in a computer 
system and a computer program product for detemiining an health indicator of a component 
at a subsequent time. A first health indicator of said component at a time, n. is detennined 
accordance with at least one coiresponding condition indicator. A second health indicator of 
the component is determined using a three state Kahnan filter at a time subsequent to time n. 



m 



In accordance with one aspect of the invention are a method executed in a computer 
system and a computer program product for ranking condition indicators used in detemiining 

10 aheal&indicatorforacomponent. Afirstsetofaplurahtyofsaidconditionindicatorsis 
detemiined. A covariance matrix corresponding to said plurality of condition mdicatois is 
detennined. A transfonnation matrix that whitens the covariance matrix is detennined. 
Differences between said first plurality of condition indicators and expected values for said 
condition indicators belonging to a health class are determined using the whitening matrix. 

15 Each health class has a conespondinghealth indicator. A portion ofsaid plurality of 
condition indicators is selected in accordance with those condition indicators have the 
smallest of said di£r»%nces. 

In accordance with one aspect of the invention are a method executed in a computer 
20 system and computer program product for estimating a conditional indicator value associated 
with a gear pair. The gear pair is modeled as a damped spring model having a contact line 
between said gears. A force. P. is detennined at a point of contact along said contact line 
causing linear and torsional response to each ofsaid two gears in said gear pair. A relative 
movement, d, is detennined of said gear pair, in acconiance with said force. P. as a sum of 
25 four responses and a contact deflection, said relative movement d representing a gear model 
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having two degrees of freedom. The relative movement, d, is used in detemuning said 
conditional indicator value for transmission error associated with said gear pair. 

In accordance with another aspect of the invention are a method executed in a 
5 computersystemandacomputerprogrampioductforestimatingaconditionindicator 
associated with a bearing. A bearing frequency ratio is determined for the bearing. A 
poiodic impulse is determined in accordance with the bearing frequency ratio. An intensity 

of an impulse on a bearing surface as a fonction of an angle relative to a bearing fault is 
detennined. A decay of a unit impulse is detennined. A movement of the bearing is 
10 determined. A conditional indicator value is determined in acconiance with the movement. 



In accordance with one aspect of the invention are a method executed in a computer 
system and computer program product for nonnalizing a set of at least one observed 
condition indicator. A plurality of conditional indicators and at least one associated factor 
15 detennined in accordance with previous data acquisitions. A mean and at least one model 
coefBcient corresponding to said at least one associated factor are determined. The set of 
observed condition indicators is adjusted in accordance with model coefficients and said at 
least one associated factor prodiicing a noxmahzed set of condition indicators. 



20 



are 



In accordance with one aspect of the invention are a method executed in a computer 
system and computer program product for determining a condition indicator about a 
characteristic of a component. A distribution of observed data associated with said 
component is detennined. A difference between said distribution and a nonnal distribution 
detennined. The condition indicator is detennined using the difference. 

25 



is 
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In accordance with another aspect of the invention are a method executed 
computer system and computer program product for detennining a condition indicator 
associated wift a component. A total impulse signal is detennined in acconlance with 
configuration data. The total impulse signal is a superposition of gear and bearing 



noise 



5 represented as a convolution of a gear and bearing signal with a gearbox transfer function. A 
condition indicator is determined in accordance with the total impulse signal. 



In accordance with yet another aspect of the invention is a method executed in a 
computer system and a computer program product for determining a health status of a 
10 componentusingatleastoneconditionindicator. At least one condition indicator is 
determined using at least one of: an impulse detennination technique and a statistical 
normality test. The health indicator is determined in accordance with the at least < 



:one 
indicator. 
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Brief Descriptioa Of Drawings 

Features and advantages of tbe present invention will become more apparent from the 
Mowing detailed description of exemplary embodiments thereof taken in conjunction with 
the accompanying drawings in which: 

5 

Figure 1 is an example of an embodiment of a system that may be used in perfonning 
vibration analysis and performing associated monitoring functions; 

Figure 2 is an example representation of a data structure that includes aircraft 
10 mechanical data; 

Figure 3 is an example of parameters that may be included in the type-specific data 
portions when the descriptor type is an indexen 

15 Figure 4 is an example of parameters that may be included in the type-specific data 

portions when the descriptor type is an accelerometer; 

Figure 5 is an example of parameters that may be included in the type-specific data 
portions when the descriptor type is a shaft; 

20 

Figure 6 is an example of parameters that may be included in the type-specific data 
portions when the descriptor type is for a gear; 

Figure 7 is an example of parameters that may be included in the type-specific data 
25 portions when the descriptor type is a planetary type; 
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Figure 8 is an example of parameters that may be included in the type-specific data 
portions when the descriptor type is bearing type; 

Figure 9 is an example of a data structure that includes analysis information; 



Figure lOisamoredetailedexampleofanembodimentofaheaderdescriptorof 
Figure 9; 

Figure 11 is an example of a descriptor that may be included in the acquisition 
descriptor groiq) of Figure 9; 

Figure 12 is an example of a descriptor that may be includc^d in the accelerometer 
group of Figure 9; 

Figure 1 3 is an example of a descriptor that may be included in the shaft descriptor 
group of Figure 9; 

Figure 14 is an example of a descriptor that may be included in the signal average 
descriptor groiq) of Figure 9; 

Figure 15 is an example of a descriptor that may be included in the envelope 
descriptor group of Figure 9; 

Figure 16 is an example of a planetary gear arrangement; 



10 



wo 02/095633 

PCT/US02/16380 



Figure 17A is an example of an embodiment of a bearing; 

Figure 17B is an example of a cut along a line of Figure 17A; 

Figure 1 8A is an example of a representation of data flow in vector transfonnalons; 

Figure 1 8B is an example of a representation of some of the CI algorithms that may 
be included in an embodiment, and some of the various inputs and outputs of each; 

Figure 19 is an example of a graphical rq)resentation of a probability distribution 
function (PDF) of observed data; 



Figure 20 is an example of a graphical rq,resentation of a cumulative distribution 
15 fimctionCCDF) observed data foUowing a gamma (5.20) distribution and the normal CDF; 

Figure 21is an example of a graphical representation of the difference between the 
two CDFs of Figure 20; 

20 Figure 22 is an example of a graphical representation of the PDF of observed data 

foUowing a Gamma (5.20) distribution and a PDF of the normal distribution; 

Figure 23 is an example of another graphical representation of the two PDFs fiom 
Figure 22 shovra which quantities as intervals rather than continuous lines; 

25 
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Figure 24A is an example of a graphical representation of the differences between the 
two PDFs of observed data and the nonnally distributed PDF; 

Figures 24B-24D are examples of a graphical data displays in comrection with a 
) healthy system; 



Figures 24E-24G are examples of graphical d^ displays in comiection with a system 
having a feult; 



Figure 25 is a flowchart of steps of one embodiment for determining health indicators 

(His); 



Figure 26 is a graphical illustration of the probabiHty of a felse almr. (f>i^^) in one 
example; 



Kgure 27 is a gr^hical iUustiation of the probability of detection (PD) in one 
example; 



Figure 28 is a graphical illustration of the relationship between PD and PFA and 
threshold values in one embodiment; 

Figure 29 is an graphical illustration of the probability of Ho and thi^old values in 
one embodiment; 

Figure 30 is an example of an embodiment of a gear model; 



10 
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Figure 31 is a graphical representation of an estimated signal having an inner beaiing 
fault; and 

Figure 32 is a graphical representation of the signal of Figure 3 1 as a fiequency, 
spectrum. 
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Referring now to Fig^ i, *own is an example of an «nbodim=.t of a system 10 
tomaybe,sedfaperf„™ngva,«ti„„ analysis and monitoring ofamactoesu* 
5 P<«<»»f»>i™*.Tl.en»ehinebeingmomtoredl2mayb.ap,rtieulardem^ 

ain:raflSensorsI4ato„^,4c are located on themactoetogato data fom one „r^^^ 
componerts of the machine. Data may be coBected by the «nsors 14. throng 14c and s«t 

Wai^ocessor or.VPUiafordatag^heringand analysis, lie VPU16analy.es ^dgathe,. 
the data from the Sensors 14a through 14c. 



The VPUI6 may also «se other data in perfcnning analysis. For example, the VITO6 

may.sec„lIec.edda.al8.0neormoreof the AlgorithmsaOmaybe used a.inp« into the 
VtVie in conn«=tion with analyzing data such as may be gafltered torn the Sensors 14. 

tough I4c. Addition.Ily.co„lig„r.,iondata22maybe.«dbytheWU16incom,ection 
with petfcmaing an analysis of the data recdved for example fiom the S««,s 14. tbon^ 
14c. Genc^y, configntati™, may include pa^metc and the lite &« may be stored in 

aeonfig„,ationdat.ffle. Each of these will be described in more detail in par.g,^hs that 

follow. 



ITK VPU16 may use as mput the con«*ed data 18, one or more of the algorithm. 20. 
and configut^ion d«a 22 to detennine one or more condition indictot. or Os. b tun, U>es. 
condition indicatorsmay be used in detennininghealth indicators or His that may be stored 

for example inaand HI storage 28. CIS describe a^tsabom.particularcomponart that 
may be usefiU in making a detemnnation about the state or health of a component as may be 
reflected m an HI depending on one or more CIs. Generally, as will be described in mo:^ 
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detail in paragraphs that follow, CIs and His may be used in connection with different 
techniques in determining an indication about monitored components such as Machine 12. 
As described in more detail elsewhere herein, the configuration data may include values for 
parameters that may vary in accordance with the type of the component being monitored. 

It should be noted that the collected data 18 may include data collected over a period 
of time fi^m sensors such as 14a through 14c mounted on Machine 12. A user, such as a Pilot 
26, may use a special service processor, such as the PPU24, connected to the Machine 12 to 
obtain different types of data such as the CI and HI values 28. 

As described in connection with Figure 1, the VPU16 may receive inputs from 
Sensors 14a through 14c. These sensors may be different types of data gathering monitoring 
equipment including, for example, high resolution accelerometers and index sensors 
(indexors) or tachometers that may be mounted on a component of Machine 12 at caretully 
selected locations throughout an aircraft. Data from these sensors maybe sampled at high 
rates, for example, up to 1 00 Mlohertz, in oixler for the VPUl 6 to produce the necessary CI 

and HI indicators. Data from these sensors and accelerometers may be acquired 
synchronously at precise intervals in measuring vibration and rotational speeds. 

Generally, the different types of data gathering equipment such as 14a-14c may be 

sensorsortachometersandaccelerometer8.Accelerometersmayprovideinstantaneous 
acceleration data along whatever axis on which they are mounted of a particular device. 
Accelerometers may be used in gathering vibration analysis data and accordingly may be 
positioned to optimally monitor vibration generated by one or more mechanical components 
such as gears, shafts, bearings or planetary systems. Each component being monitored may 
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generally be monitored using two independent sensors to provide confimation of component 
faults and to enable detection of sensor faults. 



No accelerometer is completely isolated from any other component Thus, the 
component rotational frequencies share as few common divisors as possible in order to 
maximize the effectiveness of the monitoring function being performed. For example, aU 
gears being monitored should have differing number of teeth and aD bearings should have 
differing numbers and sizes of baUs or roUers. This may allow individual components to be 
spectraUy isolated from each other to the extent that their rotational frequencies are unique. 

The indexers (index sensors) or tachometers may also be used as a particular 
monitoring component 14a through 14c to gather data about a particular component of 
Machine 12. The indexers produce a periodic analog signal whose frequency is an integer 
multiple of the instantaneous rotation frequency of the shaft that they are monitoring. These 
signals may be generated magneticaUy using one or more evenly spaced metallic protrusions 
on the shaft passing by the fixed sensor. Alternatively, these may be monitored optically 
using a piece of optically reflective material affixed to the shaft. It should be noted that each 
index point should be fixed in time as precisely as possible. In connection with magnetic 
sensors, this may be accomplished for example by interpolating the zero crossing times of 
each index pulse and similarly for optical sensors by locating either rising or falling edges. 

Assuming the minimal play or strain in the drive train when something is under load, the 
relative position and rate of any component may be calculated using a single mdex or wave 
form. 
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Because of the high data rates and lengthy processing intervals, diagnostics may be 

performed, example, on pilot command or on a predetennined flight regime or time 
interval. 



Each of the algorithms 20 produces one or more CIs described elsewhere herein in 
more detail. GeneraUy, the CI may yield useful information about the health of a monitor^ 

component TWs condition indicator or CI as weU as ffl may be used in detemiining or 
predictmg faults of different components. 



It should be noted that the VPUl 6 is intended to be used in a wide variety of 
mechanical and electrical enviromnents. As described herein, different components of an 

aircraft may be monitored. However,this is only one exampleofatype of 
which the system described hereinmay be used. .As known to those skiUed in the art. the 
general principles and techniques described herein have much broad.T and general 
appKcability beyond a specific aircraft enviromnent that may used in an example here. 



in 



m comiection with the use of CIs, the VPU16 uses the CIs as input andportions of the 
dma such as. for example, used in comrection with an algorithm to provide His. These are 
described in more detail in paragr^hs that follow. 

It should be noted that in aparticular embodiment, each mechanical part being 
monitor^ may have one or more sensors associated with it where a sensor may include for 

exampleanaccelerometeroratachometer. Generally, acceleiometers may be used, for 
example, to obtain data regarding vibrations and a tachometer may be used, for example, to 
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gain infoimation and data regarding rotation or speed of aparticular object. Data may be 
obtained and converted from the time to the frequency domain. 

A particular algorithm may provide one or more CIs. Each of the algorithms may 
> produceorbeassociatedwithaparticularCI. One or more CIs may be used in combination 
with a function to produce an HI for a particular part or type. As will be described 
detail herein, each of the algorithms may be associated or classified with aparticular part 
type. The CI generally measures vibrations and applies a function as described in acconiance 
for each algorithm. Generally, vibration is a function of the rotational frequency in the 
amom^t of torque. Using torque and a particular frequency, a CI is appropriately determined 
in accordance with a selected algorithm for apart. 

The algorithms 20 may be classified into four families or groups in accordance with 
the different types of parts. In this example, the famihes of algorithms may include shaft, 
gears, bearings, and planetary gears. Associated with each particular part being monitored 

maybeanumberofCIs.EachCImaybetheresultoroutp„tof^plyingadifferemoneof 
the algorithms for aparticular family. For example, in one embodiment, each gear may have 
an associated 27 CIs, each bearing may have 1 9 CIs, each shaft may have 22 CIs, and each 
pJanetaiy gear may have two or three CIs. It should be noted that each one of these numbers 
represents in this example a maximum nmnber of CIs that may be used or associated with a 
particular type in accordance with the nmnber of algorithms associated with a particular class 
or family. Generally, the different nmnber of CIs that may be associated with a particular 
type such as a gear by to take into account the many different ways in which a particular gear 
may fail. Thus, a CI reflects a particular aspect or characteristic about a gear with regard to 
how it may fail. 
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Different techniques used in computing CIs are described, for example, in "Introduction to 
Machinery Analysis and Monitoring, Second Edition", 1993, Penn Well Publishing Company 
of Tulsa, OK, ISBN 0-87814-401-3, and "Machinery Vibration: measurement and analysis". 
1991, McGraw-Hill Publishing, ISBN-0-07-071 936-5. 

Referring now to Figure 2, shown is an example of a data structure 50 that includes 
aircraft mechanical data. Generally, this data structure mcludes one or more descriptors 56a 
through 56n. In this embodiment there may be one descriptor for each sensor. Adescriptor 
associated with aparticular sensor includes the parameters relevant to the particular 
component being monitored. Each of the descriptors such as 56a includes three portions of 
data. The field 52 identifies aparticular type of descriptor. Each of the descriptors also 
inchides a common data portion 54 which includes those data fields common to aU descriptor 
types. Also included is a type specific data portion 56 which includes different data fields, 
for example, that may vary in accordance with the descriptor type 52. 

Descriptor types may include, for example, an indexer, an accelerx>meter, a shaft, a 
gear, a planetary gear, or a bearing descriptor type value corresponding to each of the 
different types of descriptors. The common data portion 54 may include, for example, a 
name, part number and identifier. In this example, the identifier in the common data filed 54 
may uniquely identify the component and type. 

Referring now to Figures 3 through 8, what will be described are examples of 
descriptor type specific parameters or information that may be included in a descriptor of a 
particular type, such as in area 56 of the data structure 50. 
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Referring now to Figure 3, shown is an example of parameters that may be included 
in a descriptor 60 which is an indexer descriptor type. The parameters that may be included 
are a channel 62, a type 64, a shsA identifier 66, a pulses per revolution parameter 68. a pulse 
width parameter 70, and a frequency of interest 72 for this particular type of descriptor. It 
should be noted that the type in this example for the index or descriptor may be one of 
sinusoidal, pulse such as 1/rev, or optical. The shaft identifier 66 is that as may be real or 
viewed by the indexer that calculates the shaft rate. TTxe puke width 70 is in seconds as the 
mut value. Additionally, the frequency of interest 72 for this descriptor type is a nominal 
pulse frequency that is used in computing the data quality signal to noise ratio. THe use of 
these particular data structures and parameters is described in more detail in paragraphs that 
follow. 



Referiing now to Figure 4, shown is an example of fte parameters that may be 
included in an acceletometer descriptor type 80. The descriptor for an acccleiometer type 
may include the chamiel 82. a type 84, a sensitivity 86 and a frequency of interest 88. In this 
example for the accelerometer descriptor type, the typemay be one of nonnal, or remote 
charge coupled. The frequency of interest may be used in computing the data quahty signal 
to noise ratio. The frequency of interest for a gear is the mesh rate which may be calculated 
from the gear shaft rate and tiie number of teeth of the gear. 



Referring now to Figure 5, shown is an example of descriptor type specific paramete,. 
or data that may be included when a descriptor type is the shaft descriptor. A shaft descriptor 
90 includes path parameter or data 92 and nominal RPM data 94. The path data is an even 
length sequence of gear tooth comits in the mechanical path between the shaft in question and 
a reference shaft. Tlie driving gears alternate with driven gears such that the expected 
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frequency of a gear, shaft, bearing and the like may be detennined based on an input shaft 
RPM. 



Referring now to Figure 6, shown is an example of data or parameters that may be 
5 includedinadescriptorwhenthedescriptortypeisthegeardescriptor. Included in the gear 
descriptor 100 is the shaft identifier 102 to which the gear is mounted and aparameter 104 
indicating the number of teeth in the gear. 

Referring now to Figure 7, shown is an example of an embodiment of a planetary 
10 descriptor 1 1 0 identi^g those parameters or data that may be included when the type is a 
planetary descriptor type. The planetary descriptor 1 1 0 may include an input shaft identifier 
112, an output shaft identifier 114, a parameter indicating the number of planet gears 116, a 
parameter indicating the number of teeth on the planet gear, a pararnetei J 20 indicating the 
number of teeth on the ring gear, and a parameter 122 indicating the number of teeth on the 
15 sungear. It should be noted that the number ofteeth on a planet gear relates to a planet 
carrier that is assumed to be mounted to the output shaft. Additionally, the ring gear is 
described by parameter 120 is assumed to be stationay and the sun gear 122 as related to 
parameter 122 is assumed to be mounted to the input shaft. It should be noted that the path 
between the input and the output shaft may be reduced to using a value S for the driving path 
20 tooth count and R+S as the driven path tooth count where R and S are the ring and sun tooth 

counts respectively. An example of a planetary type gear is described in more detail 
elsewhere herein. 

Referring now to Figure 8, shown is an example of a bearing descriptor 130. The 
15 bearing descriptor 1 30 may include descriptor type specific fields including a shaft identifier 
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132, a cage ratio 134, a ball spin ratio 136. an outer race ratio 1 ^8 «nH • 

, <uj ouier race ratio 138 and an mner race ratio 140 

Ai.«campteofabeari:«i,descnT«di„:„„„toail.I,ewl«r.her^ 

— and to^ ^ 
s»»Iing data a particular point in time whid, ia aaot to the VPU 16. 

ltet«ring„ow..Pig„,e9,*„™i3ancxampIe«fa,eda,ast™c,„,e.50thatcon.^ 
»«lysiaaa..Eaohi.a.anc.ofa„a,yaiada.al50aa,,pr«e«edi„*ed^ 

•headerdescriptor 152 aaddescriptorg^ups noted aal54.I„thia«amp,.tte„.„ five 
deacriptorg™^ although ^^^^ 

deaoriptor groapa 154 ttaough , iden.ifi«l ide„,fi„ ,^ j^,^ ^ 
«o«descript™.aasocia.edwia>apartioaIarg™np.^. For «a.p,e. deaeriptor g™„p ,54 
«.he acquisition groapfl^tei^.deacriptorfor each .ensorto be acquired. Tie 
acceiero^eter group >56consia.s„fadescHp,„rforeaehacce.e™ne.er.o bep,ocesaed. lie 
*«« g^up .58 includea a deacriptor for each aha« to be processed, ^e signal a«,age 
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group 160 includes a descriptor for each unique parameter set. 
includes a descriptor for each unique parameter. 

Referring now to Figure 10. shown is a more detailed example of a header descriptor 
170. Parameters that may be included in a header descriptor 170 include: an analysis 
identifier 1 72, acquisition time out parameter 174 and processing time out parameter 176. In 
this example, the acquisition, time out and processing time out parameters are in seconds. 

Refening now to Figure 11, shown is an example of a descriptor that may be included 
in the acquisition group. A descriptor 1 80 included in the acquisition group may include a 
sensor identifier 182, a sample rate parameter in Hz 184, a sample duration in seconds 186, a 
gain control setting, such as "auto" or "fixed" 1 88. an automatic gain control (AGC) 
acquisition time in seconds 190, an automatic gain control (AGQ headroom factor as a 
number of bits 192 and a DC offset compensation enable 194. 

Referring now to Figure 12, shown is an example of a descriptor 200 that may be 
included in the accelerometer group. A descriptor in the accelerometer group may include a 
parameter that is an accelerometer acquisition analysis group identifier 202, a list of 
associated planetary identifiers to be processed 204. a Kst of associated shaft analysis group 
identifiers to be processed 206, a processor identifier 208, a transient detection block size 
210, a transient detection RMS factor 212, a power spectrum decimation factor 214 specified 
as a power of 2 and a power spectrum block size also specified as a power of 2. 
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The envelope grot^ 162 
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In one embodiment, the hst of associated planetary identifiers 204 also includes two 
signal average analysis group identifiers for each planetaiy identifier, f.st identifier 
corresponding to the input shaft and a second conespondmg to an output shaft. 

5 It should be noted that the processor identifier 208 will be used in connection wift 

assigning processing to a particular DSP or digital signal processor. 

ReferringnowtoFigurel3,shownisanexampleofanembodimentofadescri^^^ 
280 that may be included in the shaft grotq,. lUe descriptor 220 may mclude a shaft ^ 

10 id«-tifi«r222.asignalaverageanalysis^upid«tifier224,ahstofgearidentifi^ 

processed226,atoofbearingidentifier.tobeprocessed228andali^ofassodate^ 
envelope analysis group identifiers 230. 

Referring now to Figure 14. shown is an example of a descriptor 2n that m.y be 
15 i-'^l-dedinthesignalaveragegroup. Itshouldbenotedthatthesignalaverage 
includesadescriptor for eachuniqueparameter set. TT.e signal averageproce 
run for each accelerometer and shaft combination even if it has the same parameters as 

another combination. Each descriptor 232 may includeanumber of output points per 
revolution 234 and a number of revolutions to average 236. 
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Referring now to Figure 1 5. shown is an example of a descriptor 240 that may 
included in the envelope g„>up. It shouldbenoted that the envelope group in^^^^^ 

descriptor for each uniqueparameter.lt is not necessary to repeat an envelopep^^^ 

each bearing if toeparameters are the same. Each descriptor 240 may includeadurati-^ 
25 parameter 242 specifying the seconds of raw data to process, an FFT size 244 which is a 
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power of 2, a lower bound frequency in Hz 246, and an upper bound frequency, also in, Hz 
248. 



Referring now to Figure 16, shown is an example of an embodiment 300 of a 
5 planetary gear arrangement. Generally, a planetaiy gear arrangement as described in , 
connection with the different types of gears and items to be monitored by the system 10 of 

Figure 1 may includeapluraUty of gears as configured, for example, in the embodiment 300. 
Included in the airangement 300 is a ring gear 302 a plurality of planet gears 304a through 
304c and of smi gear 306. Generally, the gears that are designated as planets move around 
10 the sun gear similar to that as a solar system, hence the name of planet gear versus sun gear. 
The arrangement shown in Figure 16 is a downward view representing the different types of 
gears included in an arrangement 300. 



Referring now to Figure 17A, shown is an example of an embodiment 320 of a 
15 bearing. The bearing 320 includes a ring or track having one or more sphaical or cylindrical 
elements (roning elements) 324 moving in the direction of circular rotation as indicated by 
the arrows. Different characteristics about such a structure of a bearing may be important as 
described in connection with this embodiment. One characteristic is an "inner race" which 
represents the circumference of circle 322a of the inna portion of the ring. Similarly, the 
20 "outer race" or drcmnference 322b representing the outer portion of the ring maybe a 
consideration in connection with a bearing. 

Referring now to Figure 17B, shown is an example of a cut along line 17B of Figure 
17A. Generally, this is cut dirough the ring or track within which a bearing or bearings 324 
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rotate in a circular direction. Hie ball bearings move in unison with respect to the shaft 
vsdthin a cage that follows a track as well as rotate around each of their own axis. ^ 

Referring now to Figure 18A, shown is an example of a representation 550 of 

different transfonnations that may be perfonned and the associated data flow and 
dependencies for each particular sensor. The output of the transformations are transformation 
vectors and maybe used in addition to analysis data or raw data, such as bearing ftequeacy, 
mesh fiequency, and the like, by an algorithm in producing a CI. 

Referring to the representation 550, an in going arrow represents data flow input to a 
transfonnation. For example, the FF or Fast Fourier transform takes as an input data from the 
Al signal average data transfonn. Al has as input the accelerometer data AD. It should be 
noted that other embodiments may produce different vectors and organize data inputs/outputs 
and intemiediate calculations in a variet:. of different ways as known to those skilled in flie 
art. 

Referring now to Figure 18B, shown is an example of a representation 350 relating 
algorithms, a portion of input data, such as some transformation vectors, and Cfe produced 
foreachtypeofconvonent,thatmaybeincludedinanembodiment. Other embodiments 
may use different data entities in addition to those shown in comiection with Figure 18B. As 
described elsewhere herein, each type of component in this example is one of: indexer, 

accelerometer. shaft, gear, planetary, or bearing. Certain algorithms may be used in 
comiectionwithdetemiiningoneormoreasfbrmorethanonecomponenttype. Itshould 
be noted that a variety of different algoriftms may be used and are known by one of ordinary 
skill in the art, as described elsewhere herein in more detail. The following are examples of 
some of the different techniques that may be used in producing CIs. Additionally, Figure 
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18B illustrates an example of relationships between some algorithms, aportion of their 
respective inputs and ou^uts. as well as how the algorithms may be associated with different 
component types. However, it should be noted that this illustration is not aU inclusive of all 
algorithms, aU respective inputs and outputs, and all component types. 

5 

What wiU now be described are algorithms and the one or more CIs produced that 
may be included in an embodiment. It should be noted that the number and type of 
algorithms included may vary in accordance with an embodiment. AdditionaUy. it should be 
noted that Figure 1 8B m^ not include each and every input and output for an algorithm as 

10 describedhereinandotherembodimentsofthealgorithmsdescribedgeneraUyherei^ 
also vary. 

The data quahty (DQ) algorithm 356 may be used as a quality assurance tooi for (he 
DTD CI. DQ performs an assessment ofthe raw micalibrated sensor data to insure that the 
15 entiresystemisperfoimingnominally. DQ may be used to identify, for example, bad wiring 
comiections, feulty sensors, clipping, and other typical data acquisition problems. TTie DQ 
indicator checks the output of an accelerometer for "bad data". Such "bad data" causes the SI 
to be also be "bad" and should not be used in determining health calculations. 
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What will now be described are the different indicators that may be included i 
embodiment ofthe DQ algorithm. ADC Bit Use measures the nmnber of ADC bits used in 
the cun:ent acquisition. The ADC board is typically a 16 bit processor. The log base 2 value 
of the maximum raw data bit acquired is rounded up to the next highest integer. Chamiels 
with inadequate dynamic mige typically use less than 6 bits to represent the entire dynamic 
25 range. ADC Sensor Range is the maximum range ofthe raw acquired data. This range 
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cannot exceed the operational range of the ADC board, and the threshold value of 32500 
just below the maximum pemaissible value of +32767 or -32768 when the absolute value 
taken. Dynamic Range is similar to the ADC Sensor Range, except the indicator reports 
dynamic channel range as a percent rather than a fixed bit number. Clipping indicates the 
5. number of observations of clipping in the raw data. For a specific gain vahie, the raw ADC 
bit values cannot exceed a specific calculated value. Low Frequency Slope (LowFreqSlope) 
and Low Frequency Intercept (lowFreqInt) use the first 10 points of the power spectral 
density calculated fiom the raw data and perform a simple linear regression to obtain the 
intercept and slope in the fi^iuency-amplitude domain. SNR is the signal to noise ratio 
10 observed in each specific data chamrel. A power spectral density is calculated fiom the raw 
micalibrated vibration data. For each data channel, there are known frequencies associated 
with certain components. Examples include, but are not limited to. gear mesh frequencies, 
shaft rotation rates, and indexer pulse rates. SNR measures the rise of a known tone 
(corrected for operational speed differences) above the typical minimum basehne levels in a 
15 user-defined bandwidth (generally +/- 8 bins). 

The Statistics (ST) algorithm 360 is associated with producing a pluraUty of statistical 
indicators 360a. The Root-Mean-Square (RMS) value of the raw vibration amphtude 
represents the overall energy level of the vibration. The RMS value can be used to detect 
20 major overaB changes in flie vibration level. The Peak-To-Peak value of the raw vibrating 
amplitude represents the difference between the two vibration extrema. When failures occur, 
the vibration amplitude tends to increase in both upward and downward directions and thus 
the Peak-To-Peak value increases. The Skewness coefBcient (which is the third statistical 
moment) measures the asymmetry of the probability density tunction (p.d.f.) of the raw 
25 vibration amplitude. Since it is generally beUeved that the p.d.f. is near Gaussian and has a 
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Skewness coefficient of zero, any large deviations of this value from zero may be 
indication of faults. A localized defect in a machine usually results in impulsive peaks in the 
raw vibration signal, which affects the tails of the p.d.f. of the vibmtion amplitude. The fourth 
moment (Kurtosis) of the distribution has the abihty to enhance the sensitivity of such tail 
5 changes. Khasa value of 3 (Gaussian distribution) when the machinery is healthy. Kurtosis 
values larger than 3.5 are usually an indication of localized defects. However. distributL 
defects such as wear tend to smooth the distribution and thus decrease the Kurtosis values. 

The ST algorithm may be perfonned on the following vectors: AD raw accelerometer 
10 data. Alsignal average data. RS residual data. NB narrow band data, and EV envelope data 
and others, some of which are listed in 360b. 



The Tone andBase Energy algorithm(TB) 362 uses tone energy and base en^gy. 
15 ^'^"^^"^gyiscalculatedasthesumofaUthestrongtonesintherawvibrati^ 

Localized defects tend to increase the energy levels of the strong tones. This indicator is 
designed to provide an overall indication of localized defects. "Strong tones" are determined 
by applying a threshold which is set based on the mean of all the energy contents in the 
spectrum. Any tones that are above this threshold are attributed to this indicator. Hie Base 
20 Energy measures the remaining energy level when all the strong tones are removed from the 
raw vibration spectrum. Certain failures such as wear, do not seem to affect the strong tones 
created by shaft rotation and gear mesh, the energy in the base of the spectrum could 
potentially be a powerful detection indicator for wear-related failures. Note that the sum of 
Tone Energy and Base Energy equals the overall energy level in the spectrum. 
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SI are misceUaneous shaft indicators. SOI (Shaft Order 1 in g) is the onc&-per-rev 
energy in the signal average, and is used to detect shaft imbalance. 802 (Shaft Order 2 in g) 
is the twic^-per-rev energy in the signal average, and is used to detect shaft misalignment 
GDF (Gear detector fault) may be an effective detector for distributed gear faults such as 
wear and multiple tooth cracks, and is a complement of the indicator signalAverageLl (also 
known as gearLocalFault). 

In addition to the specifically referenced vectors below, the SI algorithm takes input 
fiom the indexer zero-crossing vector (ZC). 

The Demodulation analysis (DM) 370 is designed to fiirther reveal side band 
modulation by using the Hilbert toansform on either the narrow band signal (narrow band 
demodulation) or the signal average itself (wide band demodulation) to produce the 
Amplitude Modulation (AM) and Phase Modulation (FM) signals. The procedures involved 
to obtain such signals are: 

Perform Hilbert transform on the narrow band signal (or signal average). 
Compute the amplitude of the obtained complex analytic signal to obtain the 
AM signal. 

Compute the phase angles of the analytic signal to obtain the FM signal. 
Compute the instantaneous amplitude of the analytic signal to obtain the dAM 
signal. 

Compute the instantaneous phase angles of the analytic signal to obtain the 
dFM signal. 
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The DM algorithm is perfonned on the band passed fiUered data at a frequency of interest by 
taking a Hilbert Window function of the frequency domain data and converting the data back 
to the time domain. 

The Sideband Modulation (SM) 368 analysis is designed to reveal any sideband 
activities that may be the results of certain gear faults such as eccentricity, misalignment, or 
looseness. 

CIS included in 368a are DSMn. DSMn is an indicator that characterizes the Degree of 
Sideband Modulation for the nth sideband (n = 1 , 2, and 3). The DSMn is calculated as the 
sum of both the nth high and low sideband energies around the strongest gear meshing 
harmonic. As indicated in 368b, the SM algorithm is performed on the Fast Fourier transfomi 
vector (FF). 



The Planetary Analysis (PL) 364 extracts.the AmpUtude Modulation (AM} signal 
produced by individual planet gears and compares the "uniformity" of all the modulation 
signals. 

hi general, when each planet gear orbits between the sun and the ring gears, its vibration 
modulates the vibration generated by the two gears. It is believed that when one of the planet 
gears is faulty, the ampUtude modulation of that planet gear would behave differently than 
the rest of the planet gears. The procedure to perform this algorithm is to obtain signal 
averages for the input, output, and planet shafts. For each signal average: 
Locate the strongest gear meshing harmonic. 

Bandpass filter the signal average around this frequency, with the bandwidth 

equals to twice the number of planet gears. 

Hilbert transform the bandpass filtered signal to obtain the AM signal. 
Find the maximum(MAX) and minimum(MIN) of the AM signal. 



29 



wo 02/095633 »^^„„,„-., 

PCT/US02/16380 

Calculate the Planet Gear Fault (PGF) indicator as included in 364a 
according to die equation PGF = MAX(AM) / MIN (AM). 
The inputs to the PL algorithm are the raw accelerometer data (AD) and the indexer zero- 
crossing data (ZC). 

The Zero-Crossing Indicators (ZJ) algorithm 354 is performed on the zero-crossjing 
vector (ZC). The zero crossing indicators may be determined as follows: ' 
~ -^"y+i ~ -^"y » j - - 2 , the stored zero-crossing intervals 
pulseInten>alMean=Mean{D) 

The Shaft Indicators (SI) algorithm 358 calculates miscellaneous shaft indicatora 
included in 358a. SOI (Shaft Order I in g) is the once-per-rev energy in the signal average, 
and is used to detect shaft imbalance. S02 (Shaft Order 2 in g) is the twic^per-re* energy in 
the signal average, and is used to detect shaft misalignmsnt. 

S03 (Shaft Orders), is the three-per-rev energy in the signal average, and is used to 
detect shaft misalignment. The miscellaneous shaft indicator may also be included in an 
embodiment dejSned as follows: 

p = nwnPathPairs 
p-i 

JJshqftPath^, 
. shqfiRatio:=-i=^ .^^H^ 

fishaftPath, ^""^ 



p-1 

YlmdexPath^i 
indexRatio = ^ 1- = ^^^'^^^^g 



YlindexPath,,,, 

J . n ^- indexRatio 

dnveRatio^—— pulsesUsed 

shaftRatio 



driven 
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ShaftSpeed = . ^ 



pulselntervalMean • driveRatio 
resampleRate = ^^^/^^^^ . pointsPerRev 
/?5= residual data, 

>4i = signal average, 

signaUverageLl = ^^£^ 
Rjns(Al) 

FF-FFT of the signal average, 
shaftOrderj = , j =1.3 

gearPistFauh^^^^f) 
Stdev{Al) 



As described elsewhere herein, gearDistFault (GDF) is an effective detector for 
distributed gear faults such as wear and multiple tooth cracks, and is a complement of the 
indicator signal AverageLl (also kno vwi a.>; gearLocalFaulf). 

In addition to the specificaUy referenced vectors below, the SI algorithm takes input 
fiom the indexer zero-crossing vector (ZC) and may also use others and indicated above. 

The following definitions for indicators may also be included in an embodiment in 
connection with the SI algorithm: 

shaftPath is defined for the shaft descriptor 

indexPath is the path of the shaft seen by the indexer used for signal averaging 
numPathPairs is the number of path pairs defined for shaftPath and indexPath 
pulsesUsed is the number of pulses used per revolution of the indexer shaft 
pulselntervalMean is the mean of the zero-crossing (ZC) intervals 
pointsPerRev is the number of output points per revolution in the signal average. 
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The Bearing Energy (BE) algorithm 376 perfonns an analysis to reveal the four 
bearing defect frequencies (cage, ball spin, outer race, and inner race frequracies) that usually 
5 modulate the bearing shaft frequency. As such, these four frequencies are calculated basjed on 
the measured shaft speed and bearing geometry. Alternatively, the four frequency ratios may 
be obtained from the bearing manufacturers. The energy levels associated with these four 
frequencies and their harmonics are calculated for bearing fault detection. They are: 

r Cage Energy: the total energy associated with the bearing cage defect 
10 frequency and its harmonics. Usually it is detectable only at the later stage 

of a bearing failure, but some studies show that this indicator may increase 
before the others. 

- Ball Energy: the total energy associated with the bearing ball spin defect 
frequency and its harmonics. 

1 5 - Outer Race Energy: the total energy associated with tiie bearing outer race 

defect frequency and its harmonics. 

- Inner Race Energy: the total energy associated with the bearing iimer race 
defect frequency and its harmonics. 

The Total Energy indicator gives an overall measure of the bearing defect energies. 

20 

In one embodiment, one or more algorithms may be used in determining a CI 
representing a score quantifying a difference between observed or actual test distribution data 
and a noraial probability distribution fimction (PDF) or a normal cumulative distribution 
function (CDF). These one or more algorithms may be categorized as belonging to a class of 
25 algorithms producing CIs using hypothesis tests ("hypothesis testing algorithms") that 
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provide a measure of difference in determining whether a given distribution is not nonnally 
distributed. These hypothesis testing algorithms produce a score that is used as a CI. The 
score may be described as a sum of differences between an observed or actual test 
distribution function based on observed data and a nomial PDF or nonnal CDF. An 
algorithm may exist, for example, based on each of the following tests: Chi-Squared 
Goodness of fit (CS), Kobnogorov-Smimov Goodness of fit (KS), lilliefors test of 
noimality, and Jarque-Bera test of normality (JB). Other embodiments may also include 
other algorithms based on other tests for nonnality, as known to those of ordinary skiU in the 
art. The hypothesis tests compare the test distribution to the nonnal PDF, for example as 
with CS test, or the normal CDF, for example as with the KS and Lilliefor tests. 

What will now be described is an example in which the CS test is used in determining 
a score with a test distribution of observed actual data. In this example, the test .lisrribmion 
of observed data forms a Gamma (5, 20) distribution function, having and alpha value of 5 
and a beta value of 20. The mean of this Gamma(5,20) distribution is alpha * beta having a 
variance of alpha * beta^ . The Gamma (5.20) distribution fimction is a tailed distribution 
which graphically is similar to that of a normal distribution. 

Referring now to Figure 19, shown is an example of a graphical representation 400 of 
observed data. 

Refening now to Figure 20, shown is an example of a graphical representation 410 of 
the normal CDF and the Gamma (5,20) CDF of random data. Refening now Figure 21, 
shown is an example of a graphical representation 420 of the difference between the nomial 
CDF and the Gamma (5,20) CDF. 
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In one embodiment, if there are 1 000 test samples used in fonning a single 
CDF, the graphical representation, for example, in Figure 21 represents dififerences in 1 000 
instances where the difference between the expected value (Normal CDF) and the maximum 
deviation of the (in this case defined as the score) observed gamma CDF can exceed sope 
critical value. The critical value is that statistic which represents some predefined alpha error 
(the probabiUty that the test indicates the distribution is not normal when in fact it is normal - 
this is typicaUy set at 50/0.) If the score exceeds the critical value, the distribution is said to be 
not nomial statistic. The score is the maximum deviation fiom this statistic or alpha value. 

It should be noted that the sensitivity or goodness of the test increases as the number 
of samples or instances (degrees of fi^edom V) increases apprx)ximately as the square root 
of V. For example, in the case where 1000 instances or samples ars used surdi that n-!i>DO, 
the sensitivity or ability of this CI to be used in detecting gear faults, for example, is «)ughiy 
3 1 times more powerful than kurtosis in identifying a non normal distribution. 

As another example, m the algorithm using the CS test, the normal PDF is used. 
Referringnowto Figure 22, shownis a graphical representation 430 ofthe normal PDF and 
the PDF ofthe Gamma (5,20) distribution. The representations of Figure 22 are drawn as 
continuous lines rather than discrete intervals. 



are 



Referring now to Figure 23. the quantities ofthe x-axis represented in Figure 22 
shown in anotiier representation 440 as being divided into discrete bins, intervals, or 
categories. For example, there may be 4 bms or intervals between any two integer quantities. 
Between 0 and 1. bin 1 includes values between [0,0.25), bin 2 includes values between 
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[0.25, 0.50), bin 3 includes values between [.050,0.75) and bin 4 includes values between 
[0.75, 1.0). For each bin, determine the number of observed and expected values, and their 
difference. Square each of the dififerences for each bin and then add aU the differences and 
divide by the expected value for each bin. The CS test which sums all the diJBFerences for each 
5 category divided by the expected value for each category rq>resented as: 

w ei 

for k categories or bins, k-1 degrees of freedom, fi is observed data and ei is expected data 
value or number in accordance with a normal distribution. 



10 



For each bin, take the difference between the observed and expected observation. 
Square this value and divided by expected number of observation. Sum over all bins. The 
statistic, the critical value is the at k-1 degress of freedom may be. for example, 90.72 
which is much greater than the .05 alpha value, of a x^ which is 54.57 for 39 degrees of 
15 freedom or 40 categories/bins. Thus, the observed data in this example as indicated by the 
statistic is not nonnaDy distributed. Figure 24A represents graphicaUy a difference between 
observed and expected values for each bin or interval of Figure 23. 

It should be noted that the foregoing algorithms provide a way of measuring both the 
20 skewness and kurtosis simultaneously by comparing the PDF or CDF of the test distribution 
against the PDF/CDF of a standard normal distribution in which a score is used as a CI as. 
described above. 

As known to those of ordinary skill in the art, other algorithms belonging to the 
25 hypothesis testing class may be used in computmg CIs. The particular examples, algorithms. 
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and tests selected for discussion herein are representations of those that may be included in 
the general class. 



What will now be described is another algorithm that may be used in detennining a CI 
in an embodiment of the system of Figure 1 . This may be referred to as an impulse 
detennination algorithm that produces a a indicating an amount of vibration that may be 
used in detecting a type of fault. The impulse determination algorithm takes into account the 
physical model of the system. One type of feult that this technique may be used to detect 
pit or spall on either: gear tooth, inner bearing race, outer bearing race or bearingroller 
element. This technique uses a model designed to detect this type of fault where the model 
basedonknowledgeofthephysicalsystem. For example, ifthere is a pit or spall on a 
bearing, this may produce a vibration on a first bearing which may further add vibrations to 
other components connected to or coiq)led to the bearing. 



IS a 



IS 



In one embodunent, a model caa be determined for a particular configuration by using 
configuration data, for example. In one configuration, for example, a signal received at a 
sensor may be a superposition of gear and bearing noise that may be represented as a 
convolution of gear/bearing noise and a convolution of the Gear/Bearing signal with the 
gearbox transfer fimction. Given this, ifone type offault is a pit or spall on either a: gear 
tooth, inner bearing race, outer bearing race or bearing roller element, a model that is 

designed to look for this type of fault can take advantage of knowledge of the physical 
system. 

The impulse detennination algorithm uses Linear Predictive Coding (LPC) 
techniques. As known to those skilled in the art, LPC may be characterized as an adaptive 
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type of signal processing algorithm used to deconvolute a signal into its base components. In 
the case of a pit/spaU fault, the base signal components are an impulse train generated by the 

fault hitting a surface (e.g gear tooth Avith geartooth. inner race with roUer elm^^^ 
the bearing/case transfer function. The bearing, gear and case have there own transfer 
functions. Convolution here is transitive and multiplicative. As such, LPC techniques may 
be used to estimate the total convolution function of the total vibration that may be proluced. 

For example, in this anangement. the total amomit of vibration representing the total 
impulse signal generated by a configuration may be represented as: 

in which ® represents the convoluticn operation. 

It should also be noted that convolution is a homomorphic system such that it is 
monotonicallyincreasingandthatlogariflmrictransfonnationshold. Thus the relationship of 
c = a*b also holds for Log c = Log a + Log b. A "dual nature" of convolution is used in 
following representations to equate operations using convolution in the time domain to 
equivalent multipUcation operation in the fi^uency domain. 

If "y" represents the total response of aU elementary responses, and 'V represents the 
response of the system for a series of elementary input impulses "imp" such that y is the 
convolution of imp and h, then this may be represented as: 
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and then converting "y" and "h" each, respectively, to the frequency domain 
represented as "Y" and "H", as may be represented by the foUowing: 



Y = 3(ylH = 3{h) 

taking the Fourier transform (FFT) of each where H represents the transfer function. 
The convolution in tiie time domain may be equated to a multiplication in the frequency 
domain represented as: 

Y=IMP»H 

in which IMP is the Fourier transformation of imp into the frequency domain. Above, imp i; 
in the time domain. 

The convolution in the time domain is equivalent to multiphcation in the Frequency 
Domain. Referring to the homomorphicpnq)ertyofconvolution, it follows ttiat: 

Iog(r) = log(/M>)+ log(if )^ 
therefore 

\og{lMP)=\og[Y)-\og{H\ 
/MP = exp(log(7)- log(//)) 
and finally 
imp = 3-' [imp) 

Using the foregoing, the system transfer function "H" maybe estimated for the 

Gear/Bearing and Case to recover the impulse response allocated viath a Gear or Bearing 
pit/spall fault. The estimation of this transfer function "H" may be accomphshed using 
Linear Predictive Coding (LPC) techniques. LPC assumes that the Transfer Function is a 
FIR filter, and as such, the auto-correlation of the time domain signal may be used to solve 
for the filter coefficients in a minimum sum of square error sense. 
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Using the LPC model, there is an impulse that is convoluted with a FIR fiher, such 

that: 



5 

yln] = a,x[n - 1]+ a^^" - 2]+a^x[n - 3]+ . . . 



LPC techniques maybe used to estimate the coefficients a=(aj .. an) for an order/; i 
minimum sum of square eiror sense, n=p+L. The standard least squares error estimators may 
10 be used, wherein y = y[l, 2, ...n], and x is the time delayed signal, in which: 



in a 



xin-l,n-2,.M-p] 
x[n-2,n-3yn-p-l] 



where a - (x x)" x y. These values for al .. an may be used with the following equation: 
Yhat = ax, b = (y-y/io,)^ and the estimator of error B is: 2„//b. 

15 Y may also be expressed as: 

. Y=FFr(y[l,2...n]) 

inwhich y[l.ji] aievaluesinthetimedomainexpressedinthefrequencydomainasa 
Fourier transform of the time domain values. Y represents current time vector measurements 
in the frequency domain. 

20 

In terms of a and B, the transfer function H may be estimated and represented as a/B. 
(freq. Domain). Note that "a" is a vector of the values al . . . an obtained above. 

The homomorphic property of convolution as described above may be used to 
15 estimate the impulse as represented in: 
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IMP = exp(log(Y) - log(H)) IMP Equation 

5 If there is no fault, the impulse, for example, may be characterized 

as "white noise". As the fault progresses, the impulse or the value of H becomes larger. The 
CI is the power spectral density at a bearing passing frequency for a bearing fault, or a mesh 
frequency for a gear fault Other CIs based on the foregoing value may be a "score" of the 
Lihfers test for normality, or other such test. 

10 

In the foregoing, a pit or spall may cause a vibration or tapping. Subsequently, other 
elements m contact with the ball bearing may also vibrate exhibiting behavior from this initial 
vibration. Thus, the initial vibration of the pit or spall may cause an impulse spectrum to be 
exhibited by such a component having unusual noise or vibration. 

15 

The value of IMP as may be determined using the IMP Equation above represents the 
impulse function that may be used as a "raw" value and at a given frequency and used as an 
input into an HI detennination technique. For example, the IMP at a particular frequency, 
since this the spectrum, determined above may be compared to expected values, such as may 
20 be obtained from the stored historic data and configuration data. An embodiment may also 
take the power spectrum of this raw impulse spectrum prior to being used, for example, as 
input to an HI calculation where the power spectrum is observed at frequencies of interest, 
such as the inner race frequency. For example, if the impulse ftmction is within some 
predetermined threshold amount, it may be concluded that there is no fault. 

25 

What is shown in tiie Figure 24B and Figure 24C are relative to a healthy system, 
such as a main gearbox, for example, such as in connection with a planetary race fault of an 
SH-60B U.S. Navy Helicopter built by Silorsky. 
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The Figure 24B representation 700 shows an impulse train in the frequency domain of 
the healthy system. 

5 

It should be noted that an embodiment may estimate the transfer function H usinj 5 
LPC using different techniques. An embodiment may estimate the transfer function H using 
an autocorrelation technique(AutoLPC). An embodiment may also estimate the transfer 
function H using a covariance technique (CovLPC). Use of autocorrelation may use less 
10 mathematical operations, but require more data than using the covariance. Alternatively, use 
of the covariance technique may use more mathematical operations but require less data. As 
the amount of available data increases, the autocorrelation LPC result converges to ttie 
covariance LPC result, hi one example, data samples are at lOOKKz with 64,000'data p rtlrtts 
used with the autoconelation technique due to the relatively large number of datri points. 

Figure 24C representation 710 shows the data of 700 from Figure 24B in the time 
domain rather than the frequency domain. 

20 Figure 24D representation 720 shows the power spectral density of the above figures 

as deconvolved time data of frequency v. dB values in a healthy system. 

The foregoing Figures 24B-24D represent data in a gr^hical display in connection 
with a healthy system. Following are three additional graphical displays shown in Figures 
25 24E-24G in connection with an unhealthy system, such as a starboard ring channel which 
exhibit data that may be expected in connection with a pit or spall fault. 
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Figure 24E, representation 730, illustrates an impulse train as may be associated with 
an unhealthy system in the time domain. Figure 24F, representation 740, illustrates a 
graphical display of the impulse train m the frequency domain. 

5 

In Figure 24G, shown is an illustration 740 is a graphical representation of the power 
spectrum of the impulse train represented in connection with the other two figures for the 
unhealthy system identified by a period impulse train associated with an inner race bearing 
fault. In this example, a spike may be viewed in the graphical display as well as the 
10 hannonics thereof. 



It should be noted that other algorithms and CIs in addition to those described herein 
may be used in producing CIs used in techniques in connection with His elsewhsreAerein. 

15 

What will now be described is one embodiment in which these CIs may be used. 
Referring now to Figure 25, shown is a flow chart of steps of one embodiment for 
determming the health of a part as indicated by an HI. At step 502, raw data acquisition is 
performed. This may be, for example, issuing appropriate commands causmg the VPU to 

20 perform a data acquisition. At step 504, the raw data may be adjusted, for example, in 

accordance with particular configuration information producing analysis data as output. It is 
at step 504, for example, that an embodiment may make adjustments to a raw data item 
acquired as may be related to the particular arrangement of components. At step 506, data 
transformations may be performed using the analysis data and other data, such as raw data. 

25 The output of the data transformations includes transformation output vectors. At step 508, 
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CIs are computed using the analysis data and transfonnation vector data as may be specified 
in accordance with each algorithm. At step 510, one or more CIs may be selected. Particular 
techniques that may be included in an embodiment for selectmg particular CIs is described 
elsewhere herein in more detail. At step 512, CIs may be normalized. This step is described 
5 in more detail elsewhere herein. At step 514, the selected and normalized CIs are used in 
determining His. Particular techniques for determining His are described in more detail 
elsewhere herein. 

In an embodiment, due to the lengthy processing times, for example, in executing the 
10 different algorithms described herein, HI computations may not be executed in real time. 
Rather, they may be performed, for example, when a command or request is issued, such as 
firom a pilot or at predetermined time intervals. 

The hardware and/or software included in each embodiment may vary, in one 
15 embodiment, data acquisition and/or computations may be performed by one or more digital 
- signal processors (DSPs) running at a particular clock speed, such as 40MHz, having a 
predetennined numerical precision, such as 32 bits. The processors may have access to 
shared memory. In one embodiment, sensors may be multiplexed and data may be acquired 
in groups, such as 8. Other embodhnents may vary the number m each group for data 
20 sampling. The sampling rates and durations within an acquisition group may also vary in an 
embodiment. Data may be placed in the memory accessed by flie DSPs on acquisition. In one 
embodiment, the software may be a combination of ADA95 and machine code. Processors 
may include the VPU as described herein as well as a DSP chip. 
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What will now be described are techniques for normalizing CIs in connection with 
detemining His providing more detailed processing of step 512 as described in connection 
with flowchart 500. 

5 Transmission error (T.E.) depends upon torque. Additionally, vibration depends |apon 

■j 

the frequency response of a gear. As such, the CI, which also depends upon T.E. and 
vibration, is a function (generally linear) of torque and rotor speed (which is frequency), and 
airspeed as this may change the shape of the airframe. Thus, techniques that may be used in 
connection with determining the "health state" or ffl of a component may normalize CIs to 
10 account for the foregoing since His are determined using CIs. 

For each bearing, shaft and gear wiftin a power train, a number of CIs may be 
detemiined. An embodiment may compare CI values to threshold values, apply a we?.gi itiug 
factor, and sum the weighted CIs to determine an HI value for a component at a particular 
15 time. 

Because data acquisitions may be made at different torque (e.g. power setting) values, the 
threshold values may be different for each torque value. For example, an embodiment may 
use 4 torque bands, requiring 4 threshold values and weigjits for each CL Additionally, the 
coarseness of the torque bands will result in increased, uncontrolled system variance. 
20 Alternatively, rather than use multiple threshold values and have an uncontrolled variance, an 
embodiment may use a nomialization technique which normalizes the CI for torque and rotor 
RPM (Nr), and airspeed, expressed as a percentage, for example, in which a percentage of 
100% is perfect. Use of these nonnalized CIs allows for a reduction of configuration such 
that, for example, only one threshold is used and variance may also be reduced. 
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The normalization technique that will now be described in more detail may be used in 
connection with methods of HI generation, such as the non-linear mapping method and the 
hypothesis testing method of HI generation that are also described in more detail elsewhere 
herein. 

5 

It should be noted that a deflection in a spring is linearly related to the force applied to 
the spring. The transmission may be similar in certain aspects to a large, complex spring. 
The displacement of a pinion and its corresponding Transmission Error (T.E.) is proportional 
to the torque applied. T.E. is a what causes vibration, while the intensity of the vibration is a 

10 function of the frequency response (Nr), where frequency is a function of RPM. Thus, 
vibration and the corresponding CI calculated using a data acquisition are approximately 
linearly proportional to torque, Nr.. (over the operating range of interest) and/or airq)eed 
although at times there may be a linear torque*Nr interaction effect. For example, gear box 
manufacturers may design a gearbox to have minimiun T.E. under load, and a graphical 

15 representation of T.E. vs. Torque is linear, or at least piece wise linear. It should be noted 
that test data , for example used in connection with a Bell helicopter H-1 loss of lube test, 
shows a relationship between CI and torque suggesting linearity. Additionally, tests show that 
airspeed is also relevant factor. Other embodiments may take into account any one or more 
of these factors as well as apply the techniques described herem to other factors that may be 

20 relevant in a particular embodiment or other application although in this example, the factors 
of torque, airspeed and Nr are taken into account 
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An equation representing a model minimizing the sum of square error of a measured 
CI for a given torque value in a healthy gear box is: 



C/ = Bo + B, * Torque + B^Nr + B^Airspeed + T£.. 



Equation 1) 



The order of the model may be determined by statistical significance of the 
5 coefiBcients of Equation 1 . In the previous equation, the T.E. of a "healthy" component Jnay 
have, for example, a mean of zero (0) with some expected variance. It should be noted that if 
the model fits well for the lower order. Higher order coefficients are not required and may 
actually induce error in some instances. The following example is built as a first order 
model, higher orders may be solved by extension of that explained in flie first order model. 
10 This model, written in matrix format is: y = B x where 



a, 

CI... 

a. 



B = [Bo Bu.Bn ] and x = 



1 U Nitf Airspeedl 
1 L.. Nr... Airspeed. 
1 r„ Nm Airspeedn 



Each of the CIs included in the vector y is a particular recorded value for a CI firom previous 
data acquisitions, for exanq)le, as may be stored and retrieved from the collected data 1 8. 
15 Also stored with each occurrence of a CI for a data acquisition in an embodiment may be a 
corresponding value for torque (t), Nr, and Airspeed. These values may also be stored in the 
collected data 18. 

The model coefficients for B may be estimated by minimizing the sum of square error 
20 between the measured CI and the model or estimated CI using the observed peifonnance 
data. Solving the foregoing for the unbiased estimator of B = (x^x)"' x^y . Hie variance of B 
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is: Far(B) = £(b - B)(b -B)'" = (x^x)"' where b is an unbiased estimator of B. The 

uateedestmatorofo^iss^ / =^^=fi;=£^^=^^i=5^ 

n-p-l n-p-\ n—p—l 

In the vector B from y=xB, coefficient Bo represents the mean of the data set for a 
5 particuiar component which, for example, may be represented as an offeet value. Each of the 
other values Bl . . . Bn are coeflScients multiplied by the corresponding factors, such as 
airspeed, torque, and Nr. 

The foregoing B values or coefiBcients may be determined at a time other than in real- 
time, for example, when flying a plane, and then subsequently stored, along with 
coiresponding X information, for example, in the collected data store 18. These stored values 
may be used in determining a normalized CI value for a particular observed instance of a 
Clobs in determining an HI. The normalized CI may be represented as: 

CInomatorf = T JE. = aobS - (B* X) 

where Clobs represents an instance of a CI being normalized using previously determined 
and stored B and x values. Threshold values, as may be used, for example, in HI 
determination, may be expressed in terms of multiples of the standard deviation Warning = 
Bq-^ 3*a^{x^y\ Alarm = Bo + 6*a^ (xSc)"\ It should be noted that a covariance that may be 
determined as: 

S = s^ (xSc)"^ where s^ is calculated as noted above. 

As described elsewhere herein, the foregoing techniques are based upon a healthy 
gear characterized as having noise that is stationary and Gaussian in which the noise 
approximates a normal distribution. 
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What will now be described are techniques that may be used in deteimining an HI 
using the nomialized CI values as inputs. In particular, two techniques will be described for 
detennining an HI. A first technique may be referred to as the non-linear map technique. 
5 The second technique may be referred to as the hypothesis test method of HI generation. It 
should be noted that CI values oth^ than normalized CI values may be used in connection 
with HI determination techniques described herein. 

It should be noted that an embodimmt may use CI values that are not normalized in 
10 connection with the HI determination techniques described herein. In this instance, multiple 
torque bands may be used, one for each CI or group of CIs belonging to different torque 
bands. Additionally, a largo: covariance matrix may be used as there may be a larger 
variance causing decrease in separation between classes. 

15 For any generic type of analysis (gear, bearing, or shaft), a subset of the 

diagnostics indicators or CIs is selected: The CIs which are best suited to specify the fault 
indication may be developed over time through data analysis. Faults may be calculated at the 
component level and an HI may be calculated for a given component. If there is a 
component fault, then there is a sub-assembly fault, and therefore a drive train fault. 

20 

Following is a description of a non-linear mapping methodology for determining an 
HI. Given a set of component indicators II, 12, 13, ...IN, choose the desired subset of K 
in<Kcators such that K <= N. For the chosen group of indicators, let WTi define the weight of 
the ith indicator, Wi the warning threshold, and Ai the alarm threshold. Then apply the 
25 following processing to the set of chosen indicators. 
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for XX = 1 :K /* cycle through all K indicators in subset */ 

If IPCX] < Wi /* if less than warning level Wi, assiga 0 */ 
5 Hi contribution = 0 

elseifWi*Ii<Ai 

Hi contribution = l*Wi 

else 

Hi contribution = 2*Wi 

10 end 
end 

In the foregoing pseudo-code like description, each indicator or CI is weighted and 
contributes a portion to the HI determination. Subsequently all the Hi contributions for the* 
15 selected CIs are sunomed and may be compared to threshold values for detennining one of 
two possible outcomes of "healthy" or "not healthy". 



Consider the following example table of information for a selected subset of 9 CIs 
along with threshold and weight values. It should be noted that in an embodiment, any one or 
20 more of the values for weights, warning and alarm values may be modified. 



CI 


Value 


Warning 


Alarm 


Weight 


HI contribution 


No. 




Level 


Level 






12 


3.26 


3.5 


4.0 


1.0 


0.0 


13 


3.45 


.3.0 


3.5 


1.0 


1.0 


16 


7.5 


6.0 


8.0 


1.4 


1.4 


19 


0.88 


0.5 


0.75 


0.9 


1.8 


114 


4.2 


3.5 


4.5 


1.0 


1.0 


117 


4.7 


3.5 


4.5 


0.9 


1.8 


122 


5.2 


2.0 


4.0 


1.1 


2.2 


123 


4.4 


3.5 


4.5 


1.2 


1.2 


124 


18.9 


10.0 


20.0 


1.0 


1.0 



Using the foregoing example and values, ttie sum of the HI contributions is 1 1 .4. 
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Applying the Health Indicator Contribution technique as set forth in the foregoing pseudo- 
code like description, 12, with a value of 3.26, is below the warning threshold, so the ^ 
contribution to the index is 0. Indicator 13 has a value of 3.45, which contributes a 1 toward 
the index since the weight value is also 1. However, Indicator 16 contributes a 1.4 to the 
5 index because it crosses the warning level (contributing a value of 1 to flie index) while being 
weighted by a factor of 1 .4. 

In the foregoing example, if no indicators were in alarm, the sum of HI contributions 
would be zero and if all indicators were in alarm, the sum would be 19, the worst fault case 
10 represented by this detector scheme. The HI may be represented as a value of 1 for healthy 
and 0 for not healthy as associated with a component represented by the foregoing CI values. 

The HI may be determined by dividing 1 1 .4/1 9, the maximum of worst case outcome 
to obtain 0.6. This overall health index output ratio can then be compared to another final 
15 output threshold, where normal components produce His, for example, less than 0.5; values 
between 0.5 and 0.75 represent warning levels, and values over 0.75 represent alarm. 

It should be noted that tiie weights may be determmed using a variety of different 
techniques. The wei^ts of each CI may be determined usmg any one or more of a variety of 
20 techniques. One embodiment may determine weights fortheCIsas: 

• 1 

^eigen _ values _of _ the _ cov ariance _ matrix 

It should be noted that other threshold values may be used in HI determination and 
may vary with each embodiment. 
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In one embodiment, using the normalized CI described elsewhere herein with the non- 
linear mapping technique, the threshold values may be represented as: Warning - B(? + 
3*a\xSc)'\ Alarm = Bo+ 6*a^(xSc)'\ where Bo may represent a mean or average coefficient 
5 as included in the B vector being solved for in the equations described in connection witlji CI 
normalization. In the foregoing example, the Warning threshold is 3 standard deviations and 
the Alarm level is 6 standard deviations. It should be noted that other threshold values may 
be used in and may vary in accordance with each embodiment. 

10 What will now be described is a second technique that may be used in determining 

His using CIs, in particular, using normalized CIs. 

The technique for HI deteraiination may be referred to as Hypothesis testing 
technique for HI deterannation which miniinizes the occurrence of a false alarm rate, or 

15 incorrectly diagnosing the health of a part as being included in the alarm classification when 
in fact the part is not in this particular state. In one embodiment, three classes of health 
mdication may be used, for example, normal, warning and alarm classifications with alarm 
being the least "healthy" classification. Other embodiments may use the techniques described 
herein with a different number of classes. As described elsewhere herein, the class of a part 

20 indicating the health of the part may be detemained based on measured vibrations associated 
with the part. Additionally, the technique described herein may use a transfomiation, such as 
the whitening transformation to maximize the class distributions or separation of values thus 
decreasing the likeKhood or amount of overiap between the classes. In particular, this 
maximization of class separation or distance attempts to minimize the misclassification of a 

25 part. A description of the whitening transfoimation used in herein in following paragraphs 
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may be found, for example, in "Detection, Estimation and Modulation Theory", Harry L. 
Van Trees, 1968, John Wiley & Sons , New York Library of Congress Catalog Card Number 
67-23331. 

5 Using the Hypothesis Testing method of HI generation, the HI or classification hQO 

of a vector of normalized CI values denoted as X may be detennined in which, as discussed 
elsewhere herein in more detail, X may be normalized Using the hypothesis testing 
technique, a detennination is made as to which class (normal, warning or alarm) X belongs. 
In our instance, there are three classes. However, a first determination using the hypothesis 

10 testing may be performed using a first class corresponding to normal, and a second class 
coiresponding to not normal. If the determination is normal, then testing may stop. 
Otherwise, if determination is made that the testing results are "not nonnal", a further or 
second determination using the hypothesis testing may be performed to determine v/hich "not 
nonnal" class (alarm or warning) X belongs. Thus, the hypothesis testing technique may be 

15 performed more than once in accordance with the particular number of classes of an 

embodiment. For three classes, there are two degrees of fireedom such that if the sample X is 
not firom A or B classes, then it is firom Class C. 

X may belong to class o) y or 0)2, such that: (X) {X) (the notation <> means 

20 that if q/(X) is greater than q2(X), choose class 2, o)^, or if q;(X) is less than q2(X), choose 
class 1 , 0)1 .) In the foregoing, qi is the a posteriori probability of co,- given which can be 
computed, using Bayes theorem in which qt'^ PiPi(X)/p(JO, where p(2SC) is the mixed density 
function. 
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The mixed density function is the probability fimction for all cases where q ,• is the 
unconditional probability of "i" given the probability of "i" conditioned on the mixed density 
function. 



Substituting the foregoing representation of each ql and q2, sa\CGp(X) is commoji to 

I 

o-i pAX) ^ 

both, now: P,pAX)oP^p^(X)ox as alikelihoodfimctionas W=— Trr«^- The 

likelihood ratio is a quantity in h3TJothesis test. The value P2/P/ is the threshold value, hi 
some instances, it may be easier to calculate the minus log likelihood ratio, hi this case, the 
decision rule becomes (e.g. now called the discriminate function): 

/j(A-) = -bif(A^ =-hiAW+lnp2W<>ta4- 

Assume that the piQQ's are normally distributed with mean or expected values in 
vectors Mv and covariance matrix S,- . This assumption may be detemiiaed without loss of 
generality in that, any non-nonnal distribution can be whitened, as with the whitening 
transformation described elsewhere herein, with the appropriate power transform, or by 
increasing the sample size to the point where the sample size is very large. Given this, the 
decision rule becomes: 



1^0 In ^ 



Recall that maximization of distance between the two classes is desired to mininiize 
the chance of a false alarm or misclassification of a part as broken when it is actually normal. 

A fimction Z is defined as Z = X-M, (e.g. a shift where X is the measured CI data and 
M is the mean CI values for a class), so that: 
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d^^ (z) = Z^ir^Z (this distance is the n dimensional distance between two 
distributions). 

Note that S represents the covariance. It may be determined that a particular Z maximizes the 
distance function, subject to Z^Z = /, the identity matrix. 

Using a standard Lagrange multiplier, ]x, to find the local extrema (e.g. the maximum) 
a partial derivative is obtained with respect to Z in the following: 

where S is the eigenvector of X, 

which may then be set to zero to find the extrema and solving for Z: 

, ✓ . In order that a non-nuU Z exits, Xmust be chosen to 

r''Z = fiZ or ZZ = XZ where X » '/^ 

satisfy the detenninant: ji _ a/| = o • 

Note that A is the eigenvalue of JTand S is the corresponding eigenvector. S is a 
symmetric nxn matrix (e.g. a covariance matrix), there are n real eigenvalues (X;.. .X«) and n 
real eigenvectors 0; <!)». the characteristic equation is: S4)=cDA, and <d''"<I> = I where * is an 
nxn matrix consisting of « eigenvectors and A is a diagonal matrix of eigenvalues (e.g. the 
eigenvector matrix and eigenvalue matrix, respectively). 

y, representing the coordinated shifted value of X, may be represented as: 

having a covariance matrix ofy,S^=0^2^a> = A where Sx represents the 
covariance of the vector of matrix x . Continuing, the whitening transformation may be 
defined such that: 
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Thus the transfoimation that maximizes that distance between distribution or classes 

is: 

^ =A"''^<t'^ as shown above. 
Using this value of A, define 

A'^ Si A =1, A^ S2A = K,andA^ (M2-Mi) = Land 
(S S2"^) transfoimed to a diagonal matrix A by A that may be represented as: 

A = A'^[A(I-K-^)A'']-^ A = a-K-^)'* 
which may be substituted into the discriminate function defined above: 

A(Ai=ir'-A-'r-[(-K-'x)nr+[-|i^K-'i-iHK|-in|-] 

Thus, if the above is less than the threshold, for example, ha (P2/P1), then the 
component is a member of the normal or healthy class. Otherwise, the component is 
classified as having an HI in the broken class, such as one of alarm or warning, hi the latter 
case, another iteration of the hypothesis testing technique described herein may be further 
performed to determine which "broken" classification, such as alarm or wammg in this 
instance, characterizes the health of the component under consideration. 

In the foregoing technique for hypothesis testing, values, such as the a posteriori 
probabiUties qi and qz, may be obtained and deteraiined prior to executing the hypothesis 
testmg technique on a particular set of a normalized values represented as X above: As 
known to those of ordinary skill in the art, Bayes theorem may be used in determining, for 
example, how likely a cause is given that an effect has occurred. In this example, tiie effect is 
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the particular CI normalized values and it is being determined how likely each particular 
cause, such as a normal or broken part, given the particular effects. 



It should be noted that operating characteristics of a system define the probability of a 
5 false alarm (PFA) and the probability of detection (PD). The transformation used to | 

■j 

maximize the distance function optimizes the discrimination between classes. However, the 
threshold value selected given a discriminate flmction may be used in determining the PD and 
PFA. In some embodiments, the cost of a felse alarm may be higher than the cost of a missed 
detection. In ttiese instances, the PFA may be set to define threshold values, and then accept 
10 the PD(e.g., a constant false alama rate (CFAR) type of process). The distance fimction is a 
normal density function, based on the conditional covariance of the tested vahies under 
consideration. Given that, the PFA may be determined as: Pf = P(HoHi), which means the 
probability that the sufficient statistic is greater than some threshold is the integrd of the 
threshold to infinity of a normal PDF. 

15 P« = f;'.H.(n^o)rf^=fK;,exp(-^)dr 

where 

the lower integral limit of 




In this example, the threshold may be the hi (P2/P1). This integration is the 
incomplete gamma function. Conversely, the probability of a detection (PD) is: 

but now 




25 
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Note, the distance function is relative to the condition (e.g. Uo or H/) being investigated. 

Referring now to Figure 26, shown is an example of a graphical illustration of the 
probability of a false alarm PFA represented by the shaded region A3 which designates the 
overlap between the distribution of class HO, denoted by the curve formed by line Al, and 
class HI , denoted by the curve formed by line A2. 

Referring now to Figure 27, shown is an example of a graphical illustration of the 
probability of an appropriate detection (PD) represented as area A4 as belonging to class 
represented by HI as represented by the curve formed by line A2. 

Referring now to Figure 28, shown is a graphical illusti-atlor. of a reJations^iip in one 
embodiment between the PFA and PD and the threshold value. Note that as the threshold 
increases, the PD increases, but also the PFA increases. If the performance is not acceptable, 
15 such as the PFA is too high, an alternative is to increase the dimensionality of the classijSer, 
such as by increasing the population sample size, n. Since the variance is related by 
l/sqroot(n), as n increases the variance is decreased and the normalized distance between the 
distributions will increase. This may characterize the performance of the system. The 
likelihood ratio test used herein is a signal to noise ratio such that the larger the ratio, (e.g., 
20 the larger the distance between the two distributions), the greater the system performance. 
The process of taking an orthonormal transformation may be characterized as similar to the of 
a matched filter maxunizing the signal to noise ratio. 

Referring now to Figure 29, shown is an example of a graphical illustration of how 
the threshold may vary in accordance with the probability of determining class Ho. 
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It should be noted that false alarm rate and detection rate are two factors that may 
affect selection of particular values, such as thresholds within a particular system. In the 
example embodiment described herein, false alarm rate is a deteraiining factor, for example, 
because of the high cost associated with false alarms and the fact that they may corrode 

5 confidence when a real fault is detected. It should be noted that other embodiments andj other 

"I 

applications may have different considerations. Further in this example of the system of 
Figure 1, certain factors may be considered. An acceptable false alarm rate, for example, 
such as 1 false alarm per 100 flight hours, is established. An estimate of the number of 
collection opportunities per flight hours may be determined, such as four data collections. A 
10 number of His may be selected for the system, such as approximately 800. A confidence 
level may be selected, such as that there is a 90% probabihty that a false alarm rate is less 
than 1 per 100 flight hours. 

In this example, it should be noted that each HI is a an independent classification 
event such that the law of total probability may give the system alarm rate using the 
15 foregoing: 

System PFA = 1/(100 * 4 * 800) = 3.1250 * 10 

It should also be noted that in the foregoing, when the covariance of two classes is 
approximately the same, or for example, unknown for a class, the logarithm likelihood ratio 
test for classification may be simplified in that the model may be reduced to a linear rather 
20 than quadratic problem having the following model: 
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If the covariance is whitened, the model simpUfies further (assuming the appropriate 
transfonnation is made to the means and measured values). 



What will now be described are techniques that may be used in connection with 
5 selecting a subset of CIs, such as selection of normalized CIs, for example, under 
consideration for use in determining a particular HI. 

If we have a two or more classes (such as alarm, warning and normal classifications), 
feature extraction, or determining which CIs to use in this embodiment, may become a 
problem of picking those CIs or features that maximize class separability. Note that 

10 separability is not a distance. As described elsewhere herein, an eigenvector matrix 

transfomiation may be used in maximizing the distance between two functions or distribution 
classes. However, this same technique may not be applicable when some of the information 
(e.g. dimensionality) is being reduced. For example, in the following test case, three features, 
or CIs, are available, but only two are to be selected and used in determining HI 

15 classification. The distributions are: 
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When looking at the eigenvalues of the whitening transfonnation (1.9311,3.0945, 0.4744), 
the maximum distance of the distribution is an axis y (e.g. 2"*^ dimension, the distribution was 
whitened and the project dimension (e.g. x, y or z) was plotted), but this axis has the 
20 minimum separability. Using this as one of the two features will result in higher false alami 
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rates than another feature. This may identify the importance of feature selection in 
maximizing the separabiUty. 



The problem of separabiUty may be characterized as a "mixed" problem in that 
differences in means may be normalized by different class covariance. If the mean valuejs are 
the same, or the covariance are the same, techniques such as the Bhattacharyya Distance may 
be used to measure class separability. However, same mean or covariance values may not be 
likely and thus such techniques may not be apphcable. Statistical tools developed in 
discriminant analysis may be used to estimate class separability. 

A measure of within class scatter may be represented as the*weighted average of the 
class covariance: = , for each class I, where Pi is the probability of the occurrence 

of the covariance Si for that class. In one embodiment, there may be two classes, such as 
healthy or unhealthy. When considering the unhealthy status, for example, when perfomung 
a second round of hypothesis testing described herein, there may be alarm and warning 
classes. 

A measure of between class scatta:, Sb, may be represented as the mixture of class 

means: 

Note that Mo represents the mean or expected value of the classes and Mi- Mo is a 
difference or variation from the expected value for the classes under consideration. The 
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fonnulation for a criteria for class separability may result in values that are larger when the 
between class scatter is largar, or when the within class scatter is smaller. A typical criteria 
for this is y = diagiS~'S^) , where In general, S,y is not diagonal. One technique takes the 
whitening transformation of Sh, where A^Sy4 = /, then define the whitening transformation of 
5 Sbas: 

Now taking the diagonal of the foregoing Sbw gives a better representation of flie 
class separability of each feature. 

In summary, CIs may be selected in accordance with the technique described above to 
10 obtam and examine the diagonals of tiie "whitened" Sb, represented as Sbw. l«t X be a 
matrix where rows and columns represent different CIs having a covariance matrix L. Ak 
embodiment may use normalized CIs and select a portion of these for use. An anbodiment 
may also use CIs however, those selected should belong to the same torque band. 

As described elsewhere herein, let A represent flie corresponding eigenvalue matrix 
15 and a> as the corresponding eigoivector matrix for flie CI matrix X. Then, A, as described 
elsewhere herein in connection with the whitening transformation, may be represented as: 

A = A-"2<I.^ 

where A is the transformation matrix that whitens the covariance S. If Sb is defined as above 
as the between mean covariance of the classes, the whitening matrix A may be used to 
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normalize the differences and give a distance between the mean values of the different 
classes, such that 



Sbw = A"^SbA 

where Sbw represents the "whitened" Sb. The diagonals of Sbw may then be sorted in 
descending order in which each diagonal represents an approximation of the size of the 
separation between features or CIs. Thus, selection of a subset of '*n" features or CIs from a 
possible "m" maximum CIs included in X may be determined by selecting the "n" largest 
diagonals of the matrix Sbw. In particular, the diagonal entry 1,1 corresponds to the first 
column of the covariance matrix and flie first CI in the vector X, entry 2,2 to the second 
column of the covariance matrix and the second CI in the vector X being considered, and so 
on. 

Once a particular HI is determined at a pomt in time, it may be desired to use 
techniques in connection with trending or predicting HI values of the component at future 
points in time. Techniques, such as trending, may be used in establishing, for example, when 
maintenance or rq)lacement of a component may be expected. As described elsewhere 
herein, techniques may be used in determining an HI in accordance with a vector of CI values 
having expected CI values included in vector Mi for a given HI classification, i, having a 
covariance matrix Si. One technique uses a three state Kahnan filter for predicting or 
trending future HI values. 

The Kalman filter may be used for various reasons due to the particular factors taken 
into account in the embodiment and uses described herein. It should be noted that other 
systems embodying concepts and techniques described herein may also take into account 
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other noise factors. In one embodiment, the Kalman filter may be preferred in that it provides 
for taking into account the noise of a particular arrangement of components. There may be 
noise corruption, such as indicated, for example, by tiie covariance matrices described and 
used herein. It may be desirous to filter out such known noise, such as using the Kahnan 
5 filter, providing for smoothing of data values. 

The Kalman filter provides a way to take into account other apriori knowledge of the 
system described herem. In particular, the health of a component, for example, may not 
change quickly with time. The difference between the health of a component at a time t, and 
10 time t+delta may not be large. This technique may also be used in connection with 

determinmg fiiture His of a particular part, for example, where the part is old. A part may 
have reached a particular state of relatively bad health, but still a working and functional part. 
The techniques described herein may be used with an older part, for example, as well as a 
newer part. 

15 

In the arrangement witii the Kalman filter, state reconstruction may be performed 
using the Ricatti equation, as known to those of ordinary skill in the art. The technique that 
is described herein uses a three-state Kahnan filter of HI, and the first and second derivatives 
thereof with respect to changes in time, denoted, respectively, dt^ and dt^. The Ricatti 
20 equation in this instance uses a [1x3] vector of time values rather than a smgle value, for 
example, as may be used in connection with a single state Kalman filter. 

What will now be described are equations and values that may be used in determining a 
future value of a particular HI. Let: 
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in which: 

CT is the power spectral density of the system, 

R is the measurement error, 

P is the covariance, 

g is the plant noise, 

H is the measurement matrix, 

K '\s the Kalman gain and 

O is the state transition matrix. 



H may be characterized as the Jacobian matrix. Since the value of a single HI is 
desired, only the first entry in the H vector is 1 with remaining zeroes. There are n entries in 
the n X 1 vector H for the n state Kahnan filter. Similarly, the X vector above is column 
vector of 3 HI entries in accordance with the three-state Kahnan filter. The end value being 
determined is the vector X, in this instance which represents a series of HI values, for which 
the first entry, HI_est in the vector X is the one of interest as a projected HI value being 
determined Within the vector X, represents the first derivative of HI_est and m 
represents the second derivative of HI_est. t represents the average amount of time 
between measurements or updating of the HI value. In other words, if dt represents a 
measurement or delta vahie in units of time between HI determinations, and this is performed 
for several mstances, it represents the average of the delta values representing time changes. 
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What will now be presented are equations representing the relationdiips between the 
above quantities as may be used in detennining a value of X(l) for predicting or estimating 
an HI value at a future point in time given a current HI value. 

5 

^,M='^^MM (Equation Tl) 

P„M =<I>^M|M<I>'+e (Equation T2) 

K = P^_,H^{HP^.,H^+R) (Equation T3) 

P^^iI-KH)P^., (Equation T4) 

X^=X^_,+Kim-HX^.,) (Equation T5) 

Note that the subscript notation above, for example, such as "t|t-r' refers to 
detemiiningavalueof atatimet conditioned on the measurement at a time of "t-l". 
10 Similarly, "t|t" refers to, for example, determining an estimate at a time "t" conditioned on a 
measurranoat of time "t". 

The current HI determined, for example, using otha: techniques described harein, may 
be input into Equation T5 to obtain a projected value for HI_est, the best estimate of flie 

15 current HI. To project the expect HI "n" units of time into the future, input the number of 
units of time "dt" into $ (as described above), and use the state update equation (Equation 
Tl) where now Equation Tl becomes: X t«tit = ^ X t|t . This allows the best prediction of 
HI est any number of units of time into the future where HI_est is desired. It should be noted 
that as set forth above, the linear matrix operation such as $ X is equivalent to an integration 

20 from t to.dt of the state of X, where X represents the vector of HI values set forth above. 

Different values may be selected for initial conditions in accordance with each 
embodiment. For example, an initial value for P representing the covariance may be 
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(1/mean time value between failures). An embodiment may use any one of a variety of 
different techniques to select an initial value for P. Additionally, since P converges rapidly 
to an appropriate value and the time between data acquisitions is small in comparison to the 
mean failure time, selecting a particularly good initial value for P may not be as important as 

5 other initial conditions, A value for a may be selecting in accordance with apnon 
information, such as manufacturer's data on the mean time between component parts' 
expected failure time. For example, for a transmission, the mean failure time may be 
approximately 20,000 hours. The spectral drasity may be set to (1/20,000)^. It should be 
noted that the failure rates may be generally characterized as an exponential type of 

10 distribution. The mean time between expected Mures is a rate, and the variance is that rate 
to the second power. R may also be determined using apnon infonnation, such as 
manufacturer's data, for raample, an estimated HI variance of manufacturer's data of a 
healthy component. Q may be characterized as the mean time between failures and dt (delta 
change in time between readings). As the value of dt increases, Q increases by the tibird 

15 power. 

Input data used in the foregoing trending equations maybe retrieved fix)m collected 
data, for example, as may be stored in the system of Figure 1. 

20 In determining His, for example, as in connection with the system of Figure 1 for 

particular components. His may be derived using one or more CIs. In calculating CIs, data 
acquisitions may occur by recording observed data values using sensors that monitor different 
components. There may be a need for estimating data used in connection with CI 
calculations, for example, in instances in which there may be too little or no observed 

25 empirical. For example, in connection with a power train, there may be a need to obtain 
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estimated data, for example, for each bearing, shaft and gear within the power train to 
calculate CIs. However, insufBcient empirical data may exist in connection with gear or 
bearing related measurements, such as, for example, those in connection with a gear or 

bearing related measurements, such as, for example, those in connection with a gear or 
5 bearing fault due to the rare occuirence of such events, hi such instances, mean and thre^old 

I 

values may be derived using other techniques. 

A CI may indicate a level of transmission error, for example, in which transmission 
error is a measure of the change in gear rigidity and spacing. Modehng transmission error 

10 may allow one to gauge CI sensitivity and derive threshold and mean values indicative of 
gear/bearing failure. This transmission error modeling may be referred to as dynamic 
analysis. What will now be described is a technique that may be used to model a gears to 
obtain such estimated values. By modeling each gear pair as a damped spring model with the 
contact line between the gears, transmission error may be estimated. It should be noted that 

1 5 this model uses two degrees of jfreedom or movement. Other systems may use other models 
which may be more complex having more degrees of freedom. However, for the purposes of 
estimating values, this model has proven accurate in obtaining estimates. OQier 
embodiments may use other models in estimating values for use in a system such as that of 
Figure 1. 

20 

Referring now to Figure 30, shown is an exan^le of an illustration of apair of gears 
for which a model will now be described. A force P at the contact gives linear and torsional 
response to each of the 2 gears for a total of four responses as indicated in Figure 30. The 
relative movement at P is the sum of the 4 responses together with the contact deflection 
25 due to the contact stiffoess Sc and the damping coefficient be. This may be represented as: 
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p(v 0+ '^y 2 + 

^/sp-^jcobp-mpGr /lq>'\rja}qp-Ipa> 

y + rv^/ + EQUATION Gl 

}sc+ jcobc ^ 



in which: 



sp is the linear stif&ess of the pinion; 
j is the square root of -1; 

10 (0 is the angular rate that may be obtained from the configuration file (e.g., shaft rpm 

* 60 * 2ii to obtain radians per second for the pinion driving the wheel); 
bp is the linear damping coefficient of the pinion; 
mp is the mass of the pinion; 
rp is the radius of the pinion; 

15 

kp is the angular effective stiflhess of the pinion; 
qp is the angular dampmg coefficient of the pinion; 
]^ is the angular effective mass of the pinion; 

20 sw is the linear stiffiiess of the wheel; 

bw is the linear damping coeflBcient of the wheel; 
mw is the mass of the wheel; 
rp is the radius of the pinion; 
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kw is the angular effective stiffiiess of the wheel; 
qw is the angular damping coefBcient of the wheel; 
Iw is the angular effective mass of the wheel; 

so is the linear stiffiiess of the contact patch where the two gears come into contact; 
be is the linear damping coefficient of the contact patch; 



It should be noted that values for Ihe above-refermced variables on the rights hand 
side of EQUATION Gl above, except for P (described below), may be obtained using 
10 manufacturer's specifications for a particular arrangement used in an embodiment. An 
embodiment may include quantities for the above-referenced variables in units, for example, 
such as stifi&iess in units of force/distance (e.g., newtons/meter), mass m kg units, and the 
like. 

15 The relative movement, d, is the T.E., so firom d, the above-referenced equation can 

be solved for P, the tooth force. Deflection is the force (input torque divided by the pinion 
base radius) * the elastic deflection of the shafts, which may be used in estimating P 
represented as: 

20 P = (1/kp * ip) + (1/sp) + (1/sw) + (l/kwrw)EQUATION G2 

where the variables are as described above in connection with EQUATION Gl, 

Using the above estimate for P with EQUATION Gl, the displacement, such as a vibration 

transmitted through the bearing housing and transmission case (which acts an additional 

25 transfer function), may be determined. 
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Referring again to Figure 30, shown is an example of an illustration of the gear model 
and the different variables used in connection with EQUATION Gl and G2. Lp may 
represent the longitudinal stiffiiess of the pinion and Lw may represent the longitudinal 
stif&ess of the wheel. It should be noted ttiat these elements may not be included in an 
5 embodiment using the two degrees of freedom model. 

Bearings may also be modeled to obtain estimates of fault conditions in instances 
where thei« is UtUe or no empirical data available. With bearings, a periodic impulse is of 
interest The impulse is the result of a bearing rolling over a pit or spall on the inner or outer 
10 bearing race. The intensity of the impulse on the bearing surface is a function of the angle 
relative to the fault, which may be represented as, for example, described in the Stribeck 
equation in abook by TA. Harris, 1966, Rolling Bearing Anafysis. New York: John Wiley p 
148 as: 

q(e) = qot-0^e)(l-cose)]' EQUATIONBl 

15 where n = 3/2 for ball bearings and 10/9 for rolling elements bearing, e < .5, and 9 is less than 
7c/2 in accordance wifli values specified in this particular text for the different bearings used 
in the above-refraenced Stribeck equation represented as in EQUATION Bl. 

An impulse in a solid surface has an exponential decay constant, which may be taken 
20 into account, along witii a periodic system due to rotation of the shaft The bearing model 
may then be represented as a quantity, "s", which is the multiplication of flie impulse, "imp" 
below, the impulse intensity, "q(e)" as may be determined above, the period shaft rotation, 
which is "cos(6)" below, all convoluted by the exponential decay of the material and 
represented as: 



70 



wo 02/095633 



PCT/USW/16M0. 



s^\imp^q{e)xcos0]<B>ex.v{T/t) EQUATION B2 

where T is the exponential decay and t is the time. It should be noted that 'T' varies with the 
5 material ofthesohdsmface. "exp(T/t)" may be obtained, for example, using a modal 

hammer, to generate the decay response experimentally. An embodiment may also obtain this 
value using other information as may be supplied in accordance with manufacturer's 
information. The value of "f may be a vector of times starting with the first time sample and 
extending to tite end of the simulation. T is generally small, so the expression "exp(T/t)" 
10 approaches zero rapidly even using a high sampKng rate. 

"imp" is the impulse train that may be represented as the shaft rate bearing firequncy 
ratio * sampling rate for the simulation period. 

15 "s" is the simulated signal diat may be used in determining a spectrum, "S", whore "S 

= fft(s)", the Fourier transform of s into the firequency domain fiwm the time domain. As 
described in more detail in foUowmg paragraphs, in determining a CI in connection with the 
bearing model signal "s" having spectrum "S", for example, the Power Spectral Density of S 
at a bearing passing frequency may be used as a CI. Additionally, for example, other CI 

20 values may be obtained, such as in connection with the CI algorithm comparing the spectrum 
"S" to those associated with transmission error in connection with a normal distribution using 
the PDF/ CDF CI algorithms that maybe generally described as hypothesis testmg techniques 
providing a measure of difference with regard whether the spectrum is normally distributed. 
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It should be noted that, as described elsewhere herein in connection with gear models, 
values may be used in the foregoing equations in connection with simulating various fault 
conditions and severity levels. The particular values may be determined in accordance with 
what small amount of observed data or manufacturer's data may be available. For example, 
5 in accordance with observed values, an impulse value of 0.02 for the impulse, "imp", ma^ 
correspond to a fairly severe fault condition. Values ranging from 0.001 to 0.03, for example, 
may be used to delimit the range of "imp" values used in simulations- 
Following is an example of estimated data using the foregoing equations for a bearing 
10 having the following configurations: 



Rpm = 287.1 
Roller diametCT =.25 
15 Pitch diameter = 1.4171 
Contact angle = 0 
Number of elements = 10 
hmer race fault 

20 Referring now to Figure 31, shown is an example of a graphical representation of the 

signal for the foregoing configuration when there is some type of bearing fault as estimated 
using the foregoing equations EQUATION Bl and B2. Figure 32 represents the estimated 
spectrum "S" as may be determined using EQUATION B2 above. 
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It should be noted that for bearings, there may be three types of faults, for example, 
estimated using the foregoing equations. There may be an inner race fault, an outer race fault 
or a roller element fauU. Localized bearing faults induce an excitation which can be modeled 
as an impulse train, expressed as imp in the above equation. This impulse "imp" corresponds 
to the passing of the rolling elements of the fault. Assuming a constant.inner ring rotation 
speed, the hnpulse train is periodic and the periodicity depends on the fault location. 

For outer race faults, the bearing frequency ratio, f d. or may be represented as: 



(A -/J 



EQUATION B3 



where: 

"db" represents the roller diameter, 
"dm" represents the pitch diameter, 
15 "a" =2% * frequency,yrf; 

"fir " is the rotation frequency of the inner race (e.g. shaft rate), and 
"for " the rotation frequency of the outer race (if fixed = 0). 



20 



For inner race faults, the bearing frequency ratio, f d. jr may be represented as: 



J djr 2 



r 



1+— cos(a) 



Wir-fcr) 



EQUATION B4 



Replacing a with 2ifd, the time response is f(t). This substitution may be performed 
as the initial value of a is based on an angle and not a ftinction of time. In a simulation, there 
25 is a time dependent response as expressed using f(t). 
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The radial load applied to the bearing is not constant and results in a load distribution, 
which is a ftmction of angular position. If the defect is on the outer race, the ampUtude of the 
impulse is constant because the fault location is not time varying. For an iraier race fault, the 
ann)litude with respect to angular position. The function is: 



0 elsewhere 



for \^<d^ EQUATION B5 



q(i) = q(24(0)). This quantity q(t) is amplitude at a.particular time, or q(theta) representing 
the ampUtude at a particular angle. Amplitude modulation takes into account the distance 
10 from the fault to the sensor. For outer race fenlt, the quantity cos (6) is constant (1), for inner 
race feult, it is the cosine function, noted as "cos (6)" in the above equation. 

For a linear system, the vibrations at a given frequency may be specined by the 
ampUtude and phase of the response and the time constant of the exponential decay. As the 
15 angle, 6 above, changes, the impulse response, h(t), and the transfer function H(f) also 

change due to the changing transmission path and angle of the appUed impulse. It is assumed 
that the exponential decays is independent of the angle 6, so that the response measured at a 
transducer due to an impluse appUed to the bearing at the location 9 is characterized by an 
ampUtude which is a function of 9. 

20 

The impulse response function h(t) and the transfer function H(f) may be replaced by 
a function a(9) giving the amplitude and sign of the transfer fimction H(f) 
at each angle theta and by the exponential decay of a unit impulse, (e(t)). For an inner race 
defect, rotating at the shaft frequency fs, the instantanous amplitude of the transfer function 
25 between the defect and the transducer as a fimction of time, a(t) may be obtained by 
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substituting 2it*fs*t for theta. Note that a(t) is periodic. At 0 = 0 relative to the defect and 
transducer, a(t) has its maximum value. At 6 = it, a(t) should be a minimum because the 
distance form the defect to the transducer is a minimum. Additionally, the sign is negative 
because the impulse is in the opposite direction. Because of these properties, the cos(t) may 
5 be used for the fimction a(t). 



The impulse train is exponentially decaying. The decay of a unit impulse can be 
defined by: 

10 c(0 = exp(-^^) EQUATION B6 

for f> 0, where Je is the time constant of decay. 



15 



The bearing fault model is then: 

viO = limp(tMtMO]*eit) 



EQUATION B7 



where: 



imp(t), which is the impulse over a time t, = 2ii » shaft rate * time * bearing 
frequency ratio, as may be determined using EQUATIONS B3 and B4 above; 

a(t) is the cos(e) for an inner race, which is 1 for an outer race, where cos(6) = 0, 
20 where 6 is time varying; 

and q(t) and e(t) are as described above. 



25 An embodiment may include a signal associated at the sensor for gear and bearing noise 

combined from the bearing and the gear model may be represented as: 
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sit) = [di.mtMt)a(t)]*eit)*hit) EQUATION B8 

where: 

h(t) is the frequency response of the gear case, as may be determined, for example, 
using an estimate produced with linear predictive coding (LPC) techniques or with a modal 

5 hammer analysis; 

d(t) is the signal associated with gear/shaft T.E. as may be determined using the gear 

model EQUATION Gl; 

and other variables are as described elsewhere herein. 

10 The frequency spectrum of signals representing a combined bearing and gear model 

fiom EQUATION B8 may be represented as: 

Sif) = [Dif) * F(f) * Qif) * A(J)]E(J)H{r) EQUATION 39 

As described elsewhere herein, healthy data, such as may be obtained using 
manufacturer's information, may be used in determining different values, such as those in 
connection with stiffiiesses for gear simulation, amplitude and exponential decay for bearing 
faults. In terms of geaierating fault data, since these systems are linear, the following may be 
20 defined: 

• For gear faults indicative of a crack, a reduction in the stif&ess for a tooth (e.g. 50 and 20 
percent of normal) may be used in estimating median and high fault values. Additionally, 
these values maybe varied, for example, using the Monte Carlo simulation to quantify 
variance. 

25 • For shaft misaUgnment, shaft alignment within the model may be varied to estimate 
mean fault values 
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. For gear spaliing faults, the "size" of an impulse may be detennined through trial and 

eiror, and by comparing simulation values with any limited observed fault data previously 

collected. 

• For bearing fault models, which are spalUng faults, the size of an impulse, indicative of a 
fault, with known bearing faults, may be detemiined similarly as with gear spalhng faults 

SKisitivity analysis may be performed, for example using range of different input 
values for the different parameters, to provide for increasing the effectiveness of fault 
detection techniques, for example, as described and used herein. For example, an 
embodiment may be better able to simulate a family of bearing faults to tailor a particular CI 
algorithm to be sensitive to that particular fault. 

Using the foregoing, the modulated transmission error of a gear mesh, for example, 
which is a signal may be simulated or estimated. This signal may subsequently be processed 
using any one or more of a variety of CI algorithms such that estimates for the mean and 
threshold values can then be deriyed for feult conditions. (It is assumed that the stif&iess and 
torque are known (^rion). Parameter values used in the above equations corresponding to a 
healthy gear, for example, as may be specified using manufacturer's data, may be modified to 
estimate parameter values in connection with different types of faults being simulated. By 
modifying these parameter values, different output values may be determined corresponding 
to different fault conditions. 

For example, known values for stiflfaess. masses, and the like used in EQUATION Gl 
may be varied. A cracked gear tooth may be simulated by making the stiffiiess time varying. 
The contact pitch may be varied with time in simulating a shaft aUgmnent fault. A modulated 
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inpyjX pulse on d may be used in simulating a spall on a gear tooth. Different parameter 
values may be used in connection with specifying different degrees of feult severity, such as 
alarm levels and warning levels. A particular parameter value, such as a tooth stifBiess of 
70% of the normal manufacturer's specified stif&ess, maybe used in simulating warning 
levels. A value of 20% of the normal manufacturer's specified stif&iess may be used in 
simulating alarm levels. The particular values may be determined in accordance with 
comparing calculated values with the characteristics of real CI data on any few real faults 
collected. 

While the invention has been disclosed in connection with the preferred embodiments 
shown and described in detail, various modifications and improvements thereon wUl become 
readily apparent to those skilled in the art. Accordingly, the spirit and scope of the present 
invention is to be limited only by the following claims. 
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What is Claimed is: 



1. A method executed in a computer system for determining a health indicator 

associated with a component comprising: 

determining a plurality of health classifications; 

determining at least one condition indicator quantifying a characteristic of said 
conq)onent; 

determining a probability associated with each of said health classifications, said 
probability being an estimation that said component is of a particular health classification 
given said at least one condition indicator, and 

determining, for a given set of observed values, which of said plurahty of health 
classifications is associated with said component using said probabiUties associated with said 
health classifications. 

2. The method of Claim 1, wherein said determining at least one condition indicator, 
said determining aprobabihty, and said determining which of said health clasafications is 
associated with said component, are performed upon two of said, plurahty of health 
classifications, and 

wherein, if a determination is made that a first of said health classifications is not 
associated with said component, repeating said determining at least one condition indicator, 
said determining a probability, and said determinmg which of said health classifications is 
associated with said component using two otiier health classifications. 
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5 



3. The method ofClaim 2, further comprising: 

storing previously observed and calculated data associated with said component; and 
using a portion of said previously observed and calculated data to determine said 
probability associated with said each health classification. 



4. The method of Claim 2, fiirther comprising: 

obtaining data estimates using models of said component; and 
using said data estimates to determine said probability associated with said each 
10 health classification. 

5. The method of Claim 1, fiirther con^rising: 

determining a threshold using a ratio of a portion of said probabiiilies; 
calculating a current quantity using said observed values; and 
15 comparing said current quantity to said ratio to determine which of said health 

classificatioiis is associated with said conq)onent. 



20 



25 



■80 



PCT/US02/16380 

WO 02/095633 

6. The method of Claim 5. wherein said calculating said current quantity further 
comprises: 

determining a transformation matrix maxiniizing a distance between two of said 

plurality of health classifications; 
5 determining a covariance matrix for each of said health classifications using said p 

least one condition indicator and said observed data; and 

determining said current quantity in accordance with said covariance matrices, said 

transformation matrix and said observed data. 

10 7. The method of Claim 6, wherein, for two of said pluraUty of health classifications, 

said current quantity, h(X), may be represented as: 

MAr)=iy'A-'r-[(-K-L)^]y+[-^i^K-i-|hilK|-hi|-i 

wherein, 

15 X rq)resents said at least one conditional indicator forming a vector of individual 

conditional indicator values; 
Y represents A'"^4>"^X; 

A represents (I - K'Y, diagonal matrix of eigenvalues of X using a characteristic 
equation of SO « OA, O^* = identity matrix; 
20 * represents an nxn matrix of eigenvectors, On, of X using said characteristic 

equation; 

I represents A^Si A; 
K represents A'^S2A; 
L represents A'^(M2 - Ml); 
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P2 represents an aposteriori probability of a first of said two health classifications 
given X; 

PI represents an aposteriori probability of a second of said two health classifications 
given X; 

A represents said transformation matrix maximizing a distance between distributions 
of said two health classifications represraited as A""^<&^ ; 

51 represents a covariance matrix of said first health classification having a vector of 

expected values Ml; and 

52 represents a covariance matrix of said second health classification having a vector 

of expected values M2. 

8. The method of Claim 1, wherein said at least one condition indicator indicates a 
physical state of a portion of said component. 
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9. A method for determining a health status of a component comprising: 
selecting a plurality of condition indicators having a value and each having a 

corresponding weighting factor, and at least one threshold value defining at least two 

classifications; 

determining a contribution to a health indicator for each of said condition indicatop, 
wherein said determining fiirther comprises, for each of said pluraKty of indicators: 

determining which of said at least two classifications said value of said each 

indicator belongs; and 

determining said contribution to said health indicator by said each condition 
indicator in accordance with a selected one of said at least two classifications and said 
weighted value; and 

determining said health indicator in accordance with all contributions by each of said 
condition indicator values. 



1 0. The method of Claim 9, wherein three classifications are associated with each of 
said condition indicators formed by two threshold values, and the method fijrtha- comprises: 

determining that said contribution to said health indicator for a first of said condition 
indicators is zero if said value of said condition indicator is in said first classificatioiu 

determining that said contribution to said health indicator is a first multiple of said 
20 weighting factor if said value of said condition indicator is in said second classification; and 

determming that said contribution to said health indicator is a second multiple of said 
weighting factor if said value of said condition indicator is in said third classification. 

11. The method of Claim 1 0, where the two threshold values are alarm level and 
25 warning level. 



15 
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. 1 2. A computer program product for determining a health indicator associated with a 
component comprising machine executable code for 
determining a plurality of health classifications; 

determining at least one condition indicator quantifying a characteristic of said 
component; 

determining a probabiUty associated with each of said health classifications, said 
probabiUty being an estimation that said component is of a particular health classification 
given said at least one condition indicator; and . 

determining, for a given set of observed values, which of said plurality of health 
classifications is associated with said component using said probabilities associated with said 
health classifications. 

13. The computer program product of Claim 12, wherein machine executable code 
for determining at least one condition indicator, determining a probabiUty, and determining 
which of said health classifications is associated with said component, are perfonned upon 
two of said plurality of health classifications, and the computer program product fiuAer 
mcludes machine executable code, wherein, if a determination is made that a first of said 
health classifications is not associated with said component, for repeatedly determining at 
least one condition indicator, repeatedly determining aprobabiUty, and repeatedly 
determining which of said health classifications is associated with said component using two 

othor health classifications. 

14. The computer pro-am product of Claim 13, further comprising machine 

executable code for: 

storing previously observed and calculated data associated with said component; and 
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using a portion of said previously observed and calculated data to determine said 
probability associated with said each health classification. 

15. The computer program product of Claim 13, further comprising machine 

5 executable code for: 

obtaining data estimates using models of said componart; and 
using said data estimates to determine said probabiUty associated with said each 
health classification. 



10 



16. The computer program product of Claim 12, further comprising machine 

executable code for: 

determining a threshold using a ratio of a portion of said probabilities; 
calculating a current quantity using said observed values; and 
comparing said current quantity to said ratio to deteraiine which of said health 
15 classifications is associated with said component. 

17. The computer program product of Claim 16, wherein said machine executable 
code for calculating said current quantity fiurther comprises machine executable code for. 

detennining a transformation matrix maximizing a distance between two of said 
20 plurality of health classifications; 

determining a covariance matrix for each of said health classifications usmg said at 
least one condition indicator and said observed data; and 

detennimng said current quantity in accordance with said covariance matrices, said 
transformation matrix and said observed data. 

25 
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18. The computer program product of Claim 17, wherein, for two of said plurality of 
health classifications, said current quantity, hQQ, may be represented as: 

M;0=|i''A-'y-[(-K-'i)^F+[— I^K-'l-|hi|K|-hi|.] 

5 wherein, 

X represents said at least one conditional indicator forming a vector of individual 

conditional indicator values; 
Y represents A*''^«>^X; 

A represents G - K*') diagonal matrix of eigenvalues of X using a characteristic 
10 equationofS<I> = <I>A,«'^'I'= identity matrix; 

* represents an nxn matrix of eigenvectors, «I>1.. «I»n, of X using said characteristic 

equation; 

I represent^ A^BiA; 
K represents A^S2A; 
15 L represents A''^(M2- Ml); 

P2 represents an aposteriori probability of a first of said two health classificatioiK 

givenX; 

PI represents an aposteriori probabiUty of a second of said two heallh classifications 
given X; 

20 A represents said transformation matrix maxunizing a distance between distributions 

of said two health classifications represented as 
A-'V; 

Si represents a covariance matrix of said first health classification having a vector of 
expected values Ml ; and 
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S2 represents a covariance matrix of said second health classification having a vector 
of ecpected values M2. 

19. The computer program product of Claim 12, ydierein said at least one condition 
indicator indicates a physical state of a portion of said component. 
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20. A computer program product for determining a health status of a component 

comprising machine executable code for: 

selecting a plurality of condition indicators having a value and each having a 

corresponding weighting factor, and at least one threshold value defining at least two 
5 clasafications; 

determining a contribution to a health indicator for each of said condition indicators, 
wherein said machine executable code for determining further comprises, for each of said 
plurality of indicators, machine executable code for: 

detemuning which of said at least two classifications said value of said each 

10 indicator belongs; and 

determining said contribution to said healtii indicator by said each condition 
indicator in accordance with a selected one of said at least two classifications and said 
weighted value; and 

determining said health indicator in accordance witii aU contributions by each of said 
15 conditioQ indicator values, 

21. The computer program product of Claim 20, wherem three classifications are 
associated with each of said condition indicators formed by two threshold values, and tiie 
computer program product further comprises machine executable code for: 
20 determining that said contribution to said health indicator for a first of said condition 

indicators is zero if said value of said condition mdicator is in said first classification; 

detennining that said contribution to said health indicator is a first multiple of said 
weighting factor if said value of said condition indicator is in said second classification; and 
determining that said contribution to said healtti indicator is a second multiple of said 
25 weighting factor if said value of said condition indicator is in said fliird classification. 
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22. The computer program product of Claim 21 , where the two threshold values 
alarai level and warning level. 
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23. A method executed in a computer system for determining an health indicator of a 
component at a subsequent time comprising: 

determining a first health indicator of said component at a time, n, in accordance with 
at least one corresponding condition indicator; and 

using a three state Kahnan filter to determine a second health indicator of said 
component at a time subsequent to time n. 

24. The method ofClaim 23, further comprising: 

estimating said first health indicator using a hypothesis determination technique. 

25. The method of Claim 24, wherein the health of a component does not vary by a 
large amoimt with respect to a change in time. 

26. The method of Claim 24, wherein said at least one condition indicator is a 
normahzed condition indicator. 
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27. The method of Claim 23, wherein 
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// = [l 0 o] 



1 dt ^% 



<D= 0 1 dt 
0 0 1 




Eldest 
HI 



HI 



5 in which: 



a is a power spectral density, 

R is a measurement error, 

P is a covariance, 

g is a plant noise, 

H is a measurement matrix, 

ATisaKahnangainand 

O is state transition matrix, and 



and the method comprising: 

using the above Equations Tl through T5 to determine a value Hi_est which is a 
estimate of a value ffl at a time in the future. 



(Equation Tl) 
(Equation T2) 
(Equation 13) 
(Equaction T4) 
(Equation T5) 
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28. A computer program product for determining an health indicator of a component 
at a subsequent time comprising machine executable code for 

determining a first health indicator of said component at a time, n, in accordance with 
at least one corresponding condition indicator, and 
5 using a three, state Kahnan filter to determine a second health indicator of said 

component at a time subsequent to tinie n. 

29. The computer program product of Claim 28, further comprising machine 
executable code for: 

10 estimating said first health indicator using a hypothesis determination technique. 

30. The computer program product of Claim 29, wherein the healA of a component 
does not vary by a large amount with respect to a change in time. 

15 31. The computer program product of Claim 29, wherein said at least one condition 

indicator is a normalized condition indicator. 



20 



25 



92 



wo 02/095633 



PCTAJS02/16380 



32. The computer program product of Claim 28, wherein 



10 



15 



0 = 



in which: 

cy is a power spectral density, 

R is a measurement error, 

P is a covariance, 

Q is a plant noise, 

jFf is a measurement matrix, 

a: is a Kahnan gain and 

<D is state transition matrix, and 



[1 0 


0) 






1 dt ^^y^ 

0 1 dt 




0 0 


1 






dt 1 








m 




jest 






m 






m ■ 







X^=X^_,+KiHI-HX^.,) 



(Equation Tl) 
(Equation T2) 
(Equation T3) 
(Equation T4) 
(Equation T5) 



and the computer program product comprises machine executable code for: 
determining, using above Equations Tl through T5, to determine a value Hi_est 
20 which is a best estimate of a value HI at a time in the future. 
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33. ^"'^-^^^^^^^f^ 

^ '^l' of a«i differences havim , 

"»«^iida«condiiion indicator 



con»pomltoanobs.,v«idaU.«p^B„^ 



accordance v/ith those 
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36. The n>clh„d Of CW™ 33. ftote c«priAg: 



part.c„.arcIass.aodK«fl»p„*^„„„,.„,^. 

37. »«'n«Wofaaim36,firt«r«m,p™i„g: 
S^-A'Syi. for ae «hiteDi,« traisfemMion mate A. 
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WO02/095«3 

A ^„ f CTAJS02/16380 

ahealthindicatorfcracn™. . "«<1 u> drtamofflg 

„, ,,, ^ ^"^^''"^wnipriang machine 
racecutable code for 

»2^----'»<.c^o^.ea.o....«^,^^/ 
coiresponding condition indicator. 



42. Th^owwuterprogrm product of Claim 39 fimh. 
—.ecodefe ^'-^'•«^«'»pn™«m«*i«e 

de«« a mea^e Of 5e,weeo Cass ^. S.. .ep.^.^ 
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^ , . '-'*^''*«'««nte»expectedval„eofa 

> executable code for: 

^^--ii^. for tte white^g ^fo^^ ^ 

fransfonnationinatrixAis; ^ 
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comprising: 

gears; i 
■"""■^ «<* "f^a two 6«r3 fa said gear pai; 

. gear modd having two degrees of fi«dom; and 
•rawmssion em>r aasodaled with said gear pair. 



46. Tie "■«lK>dofClaim45,fi«hercopvriai»g: 

v»>^p.rao«erv.aeai„aee„rd..e.withsi™.,ati.gditr^^ 
d«e™i«i„g„ea„3odth^„„^^,„^,^^,^^^^^ 
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47. A method executed i„ a computer system for estimating a conditional indicator 
value associated with a bearing comprising: 

determining a bearing frequency ratio for said bearing; 
determining a periodic impulse in acconiance with said bearing frequency i^o; 
5 determining an mtensity of an impulse on a bearing surface as a fimction of an angle 

relative to a bearing fault; 

detemiining a decay of a unit impulse; and 
determining a movement of said bearing; 

detennining a conditional indicator value associated with said bearing in accordance 
> with said movement. 



48. The method of Claim 47. wherein said bearing frequency ratio corresponds to 
of: an inner race frequency and an outer bearing race frequency. 



one 



49. The method of Qaim 47, fimher cornprising: 

vaiyingparameter values in determining conditional indicator vahies associated with 
different health classifications of said bearing. 

50. Hie method of qaim 49, further conq>rising: 
estimating values associated with a gear noise; 
estimating values associated with a bearing noise; and 

estimating a condition indicator of a gear and bearing noise as a combined signal 
using said values associated with said gear noise and said bearing noise. 

51 . The method of Claim 50, fiirther comprising: 
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determining said combined signal s(t) as: 

« [rf(0/(/)9(0fl(0]* e(0 * A(0 
where: • 

5 is a frequency response ofa gear case; 

^r?^ is a signal associated with gear and shaft transmission enor, 
fl[t) is a bearing frequency, ratio; 

q(t) is an amplitude at aparticular time t, 

.(.)fe.cosi„eibroncofantaerand™ter,.ceconditoa,a,<,c«i„„«^«^^ 
an impulse is applied to said bearing; and 

e(t) is a decay rate of unit impulse; 

52. -me meUKK, of data 51, whctin h(,) is toe^ined using c»« „£ a linear 
pi«B<*ye coding (IfQ tecWqoe and moM lammsr aiBlysis. 

53. ThemMhodofClaimSI.ftothercomprising: 

'tete-Mneafie,»„ey,,«ot™ofsig^ 
by forming: 



m the frequency domain. 



54. A computer progr^ product for estimating a conditional indicator value 
associated with a gear pair comprising machine executable code for: 
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gears. 



and 

gearpain 

determining a relative movement, d. of said i 



Wona] response to each ofsaidtwogea.i.sa,dgearpai. 



gearpair, in accordance with saidforce 

using said relative movement d in •• 

^'"^''"^'^^"^Ssaidconditionalindicatorvaluefor 
transmission error associated withsaidgearpair. 
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IS 



55. The computer programproduct of Claim 54 forther . 

v.jami 54, fiirther compnsmg machine 

executable code for: 

relative to a bearing fault; lanangie 
detemiining a decay of a um't impulse; and 
determining a movement of said bearing; 

deteimimngaconditionalindicatorvalueassociatedwithsaidb • 
withsaidmovement. beanng m accordance 
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57. The computer program product of QaimSfiu,!, 
con^ndstooneofaninn ^'-^—<'''«anng frequency nuic 

oneof.anu.erracefrequencya„d.outerbe^g,.eeiiequ.^^^ 

58. The computer program product of Claim ^fifi^i, I , 

executable code for: f 

executable code for: 

vdues «soci.ted . gear 
««n,att,g values ^«i,ted wia a b«ri„g ^ 

'^'"^''»'^-'0^.*-»d^..„ia.a.a«^^,,, 

60. ThecomputerprogramproductofCIaim59 fi^h 
executablecodefon comprising machine 

determining said combined signal s(t) as: 



where: 

*ry is a fi«,„ency r=^^ of a gear ca«; 

is a si^ a^ociated gear and shaft ^^^^ 
f(t). is a bearing frequency ratio; 
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q(t) is an amplitude at a particular time t. 
a(t) is a cosine for one of an 
an impulse is applied to said bearing and 
eCt) is a decay rate of unit impulse; 



innerandouterraceconditionatalocationtheta at Which 



- ^^-P--.amp.ductofCla^.o,.here.h(t)isd^^ 
of:ahnearp.edictivecodingtechm,ue.d modal hammer analy.3. 



using one 



represented as 



m 
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63. Amethod nonnalizing a set of at least one observed condition 
comprising: 



indicator 



5 "caWance with previous data acqmshioiis; 

de,e™ni„gamea„andatleastone„,«idcoefBci«rtam«po„di^ 
one associated factor, and 

adjusdng said se. of *s«ved coKUti™ i,tf oata« i. accori^c. 
oocffidcts and said a. leas, on. associaW fiCor pn^udng a n«^ 3e, „f condition 

10 indicators. 



64. The method of Claim 63, wherein said relevant fecto« includ. at k^t 
torque, airspeed, and rotor speed. 



one oi- 



ls 



65. ^-ethodofClaim64,whe.insaidrelevantfactorsincludeton,ue,a^ 
and rotor speed (Nr), and the method further comprising: 

solving forB using the equation form y=B X for the following: 



y= 



c/, 

A. 



B = [ffo Bu..Bn ] andx^ 



1 ti iVs, Airspeedl 
1 Nr... Airspeed... 

} t„ Nm Airspeedn 



whetdn .i>e vector y includes condidon indicate values and the matrix x indudes relevant 
20 ^"^"'^"""^"fing to one ofsaid condition indicators included in y. and Binelnte 



coefficients. 
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66. The method of Claim 65 wherein o 
bme by deteimining: ^ 

CInonnahzed = dobs - (B * x) 



.0 """^"^""•^"^"^——e ™^SC.„™^ 



oint 



'^'^ "P-^ .-hi.. ™. ^ ^ 



69. ThecomputerprogramproductofClaimfiR^i, 
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uy,ue. airbed a.dn>,or.peedW.a^fl„eo™pu,„^^^^^^^ 
machine executabJe code for: 

solving for B using the equation f(mn y= B X for the foUowing: 

} K, Airspeedn 

coefficients. 



71. •n"='»Wuterprogrmp,od«ctofaaim70,,>d«a.,^^^CIi, 
dc.cn™^a.a.,.crpoi«i„«.e.a„dfl«c<»,„,erTOp™,«fi^„^^ 
machine executable code for detennining: 

Qnoimalized = Qobs - (B * x) 

wh<«CIobsisana<«<*sen,«ivaI„ccon«p«,di.^to.,K,^^ 

Clhonnalized. 



72. neoon^programproduaofClaim7I.iteh«co„^^^ 

executable code for: 

storing said values for B and x in a database; and 

retrieving said values forBandxinconnection with detennining 
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chararteristic of a component comprising: 

detent a di«»«i™ of ob^rved da« a^aW 

74. The melliod of CWm 73, fimher comprising: 
de.«^g^ete«Udd,«„^„,oaofofcse™=dda^i.^^^ 

noimaKty. 



with said at least one normality test. 

76. ^ of Claim 75. „.^b««oo of obsen-ed daU a«ociated 

«a a compone.. app^xima^ oae of: a Gaussian d.-3^^^^^ 
and a non-Gaussian distribution otherwise. 
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77. The ir^ethod of Claim 74. iurth^- comprising: 

sa.dexpecteddatabeing.epresentedbysaidnonna,di.tribution;^ 
determining a sum using said differences; and 

^^-----.-e^cesisgreatertbanacri^valu, 
observeddataisnot„onnaIWs,ribute.saidcriticaIvaluebeinga 

78. The method of Claim 77. further comprising: 

detenniningascorebeingamaximumde.-ahonfK».s^^ 
condition indicator being said score. 



a» a nmnlw of observed data values iacnases. 



80.^„e«.odofCIai,.,3.«^ei„^^^^^^,^_^ 
distribution of expected values. 
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81. A computer program product for detennining a condition indicator about a 
characteristic of a component comprising machine executable code for: 

detennining a distribution of observed data associated with said compoiient; 
measuring a difference between said distribution and a nomial distribution; Jpd 
detennining said condition indicator using said diflference. 



82. The computer program product of Claim 81, further comprising: 

machine executable code for detennining whether said distribution of observ^^ 
normally distributed using said difference using at least nomiality test that is one of: chi- 
square goodness of fit test, Kohnogorov^mimof goodness of fit test, Lilliefois test of 
normahty and Jarque-Bera test of nonnality. 

83. The computer program product of Claim 82, wherm said nonnal distribution is 
one of a nonnal cmnulative distribution function and a normal prx,bability distribution 
function in accordance with said at least one nonnality test. 

84. The computer program product of Claim 83, wherein said distribution of 
ob^eddataassociatedwitfaa component approximates one of: a Gaussian distribution if 
said component is healthy and a non-Gaussian distiibuh-on otherwise. 



110 



10 



wo 02/095633 

PCTAJS02/16380 

85. 'n»ecomputerprogrampnKluctofClaim82,lurtherconiprismg^^^^ 
executable code for: 

determining a number of differences between said observed data and expected data, 
said expected data being represented by said normal distribution; and 
5 determining a sum using said differences; and 

determining that said observed data is not normaUy distributed, said critical val« 
being a threshold if said number of differences is greater than a critical value. 



86. The computerprogram product of Claim 85, further comprising: 
machine executable code for detemiinmg a score being a maximum deviation fiom 
said critical value, said condition indicator being said score. 



87. The computer program product of Claim 86, wherein sensitivity of said condition 
indicator increases as a number of observed data values increases, 

15 

88. The computer program product of Claim 81. wherein said normal distribution 
^proximates a distribution of expected values. 
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89. Am«hod«<«di,aco™prtersystemf„rd=ta»iiBgaco^^^ 
asKraated wilh a component, Ihe method comprising: 

detemmung.total impulse «g™, in acc«da„^wfflo«tfg„n^„,d^^,to^ 

imp»lsesig,«lbeingasup.,posi.ionofgeara„dbe»ng„oiser,™,s.con™^ 
5 ofageKtmJbeaiingsignalmtii.geaitoxtransfertoctionia^ 

aelemdning a cozrfition indicator in accontance »ith said total itW^^ 

90. Themeaod of Qaim 89, further comprising: 
n,»«»«ng a total in^ sig^ by . 

said component as: 
in Which ® represents a convolution operation. 



15 



20 



91. The method of Claim 90, fortha- comprising: 

representing convolution operations inatime domain to equivalent op^^ 
frequency domain. 



9Z The method of Claim 90, finther comprisiiig: 

asatransferftnctioninafi^p^ 

domain nsing a linear predictive coding techm-qne to deconvolve a signal into it. base 
components. 



25 

93. The methbd of Claim 92, further comprising: 
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estimatingsaidtransferfu„ction,H. in said frequency domain as a/B, wherein a = 
(al, ... an), each ai representing an ith coefficient for an order p. n=p+l. as: 

yl»ha,xln-l]+a,xln-2]+a,x[n-3]+... 



and B is an estimate of an eiror represented as: 
in which: 

10 »> = (y-yU', y=y[l,2,...n], 

yAaf is an estimated value of y, yj^ = a:^ 
X is a time dekyed signal rqiresented as: 

^x[n ~l,n—2,.ji-p] 
Jt[n-2,»-3,«-p-lJ 



15 x= 



where a = (x'" x) ' x^ y, values for al 



an. 



20 



94. The method of Claim 93, further compiisii^g; 

estimating an impulse, D^. to said fi«juency dom^ 

IMP = exp(log(Y) - log(H)), 
in which: 

Y«ffl(y)andH»ffl(h),wheremistheFouriertr^fon„fi„,ction,y^^ 
time domain, Y and H are in said frequency domain. 



25 



95. ThemethodofClaim94,whereinavalueassociatedwithH 
increases. 



mcreases as a &ult 
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96. Tl.en.ethodofClaim94,wh9«„stid«mditi»j„dicatorisaaid,alueofD^^ 



W. TlK mahod of aaim 94, itofta comprising: 

condition indicator being said vaJue. 

98. ThemethodofClaimPT.whereinsaidfrequencyofinterestisatl^ 
bearingpasskgfi^quencyforabearingfault,andameshfi^uencyforag 

99. The method of Claim 98, further comprising: 

perfonningaFouriertransfonnationtoobtainllVJPinsaidfi^i^^^^^^^^ ^ 

100. The method ofClaim 89, further comprising: 

detectingafaultinconnectionwithpredetenninedvaluesofsaidhealthstatus^^ 
said conditionindicator,v*erein said faultbeingdeteotedis«,eof^ 

ageartooth,innerbeariBgxace,outerbearingrace,andbearingroUerelement 
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detenmoiog a U«al tapulsc signal in acoortanoe wi^ 

of a gear aiid bearing sisiBl with a geaAox taiBfe ftnctioii; and 
deh^ining a «™ii«„n indio*,r in accoriance «tt 

IM. Thecon^pr„g^p^„fClaim lOI.flmhcrcompriaingmac^ 

executable code for: 

«prc«nfagaWalin,p„,aeagnalgene™edbyaconfig„,aUonofaa^ 

said component as: 

in which erepresenlsa convolution qxraSon. 



103. •n^con.p^erprogram ofClain.102. <u«her comprising „^„^,, 

code for: 

^^P^sentingc^^Monopemionainatimedon^to.^^^^ 
frequency domain. 
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104. 'I^^^^o^Pmerprogram product of Claim 102, furfter^^^^^ 
executable code for: 

5 components. 



105. The conq,uterprogram product of Claim 104, further comprismg madiine 
executable code for: 

estimating said transferfunction,H, in said frequency domain as a/B. wherein a = 
10 (al. .., an), eadi ai rq)reseating an ith coefficient for an order p,n=p+l, 

){nha,xln-l]+a^jin~2]+a,xln~3]+... 



as: 



and B is an estimate of an error represented as: 

in which: 

»>=(y-y&«)^ y=y[i,2,...n], 

20 Yhat is an estimated value of y, y^, = aj^ 
X is a time delayed signal rqjresented as: 



x= 



Jc[«-2,n-3,«-/>-l] 



where a = (x^ x)"' x^ y, values for al ..an. 

25 

106. ThecomputerprogramproductofClaim 105, further comprising machine 
executable code for: 
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component as: 



10 



estimating an impulse, IMP, in said foquency domain of said 

IMP = exp(log(Y)-Iog(H)), 
in which: 

Y - ffl(y) ai.d H = ffl(h). where fft is the Fourier transfonn function, y ajid h are in a 
time domain. Y and H are in said frequency domain. 

107. The computerpiogram product of Qaim 106, wherein a value associated with H 
increases as a fault increases. 

108. ThecompmerprogiamproductofCIaim 106, whereinsaidconditionin^^^^^ 
is said value of IMP. 



109. The computer program product of Claim 106, further comprising machine 
15 executable code for: 

calculatbg a power s^al <ta^,y of said in^ulse JMP in a fiequ^ 
detennimng a val« of said power apecW density at a flequeocy of interest, said 
condition indicator being said value. 



20 ""•^«»^»'P"8>™I»^ctofaaiml09,whe,™isaidfte,uencyo 

is a. least one of: a bearing passing fl«,„«cy for a bearing fault, and a n»sh teq»^ 
gearfault 



111. The computerprogram product of Qaim 110, further comprising machine 
25 executable code for 

perfomiing a Fourier transformation to obtain IMP in said frequency domain. 
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112. Thecomputea-programproductofClaimlll.lurthercomprisingmac^^ 
executable code for; 

detecting a fault in connection with predetennined values of said liealth status using 
said condition indicator, wherein saidfault being detectedisoneofapit and s^^ 
a gear tooth, inner bearing race, outer bearmg mce. and bearing roUer el^^^^ 
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113. Amethodexecutedinaccanpatersystemfordetenniningahealth^^^ 
component using at least one condition indicator, the method comprising: 
detennining said at least one condition indicator using at least one of: an impulse 
determination technique and a statistical normality test; and 

detennining said health indicator in accordance with said at least one condition 
indicator. 



114. ra=n>*odofci™il3,whmtosaidstatistiod«o™»n,yh«fao»eol!chi. 
sq»»e goodn«s of fit ,«t, Kotoog,«^^,^f of fit tea. lilliefor. test of 

normality and Jarque-Bera test of normaUty. 

115. •n>'«'=tt»dofCIaimU3.„hereinexp«=tertdaUvrf««w™ 
distribution. 
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116. Acomputerprogramproductfordeterrniningahealthstatusofa^^^^^ 
using at least one condition indicator, the method comprising: 

detenniningsaidatleastoneconditionindicatorusingatleastoneof: animpulse 
detennination technique and a statistical normality test; arid 

detemuning said health indicator in accorfance with said at least one JditiJn 
indicator. 

117. The computerprogram product of Claim 116, wherein said statistical nonnality 
test is one of: chi-square goodness of fit test, Kolmogorov-Smimof goodness of fit test. 
Lilliefors test of nonnality and Jarque-Bera test of noimality. 

118. TTiecomputerprogramproductofCiaimll6,whereine,q>ectedd«tavalues 
approximate a normal distribution. 
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